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Abstract: Superstring loops formed by intercommutation of low tension horizon- 
crossing superstrings may be captured and accreted by growing matter perturbations. 
The paper explores the influence of string tension fi and of the network formation mecha- 
nism on the clustering process. Galaxy formation and growth is schematically described 
by a radial infall model. A fully relativistic treatment of motion in perturbed Friedmann- 
Robert son- Walker cosmology is developed and applied to track the loop center of mass 
motion. The local enhancement of loops (galaxy's loop density divided by universe's 
loop density) is calculated and compared to that of cold dark matter for different fi in 
the context of two network formation scenarios, one yielding large loops (fragmentation 
near the horizon scale) and the other small loops (cusp- mediated formation). 
The primary physical process that enables capture of a loop by a growing perturbation 
is the decrease in peculiar velocity slaved to the universe's expansion. The main process 
that removes the loop from the galactic potential well is the "rocket effect", the recoil 
of anisotropic gravitational wave emission. Quantitative criteria for capture and detach- 
ment are given. There is a critical value of /i below which clustering generically occurs 
in our own Galaxy. 

Fragmentation scenarios for large-scale loops lead to significant halo enhancements. With 
typical model parameters (velocity dispersion of newly formed loops, loop length distri- 
bution, etc.) the limiting enhancement of loop (energy) density is ~ 0.25 — 0.4 that of 
cold dark matter which is fully achieved for Gfi/c^ < 10~^^ at all scales ^10^ kpc. The 
fact that the string loop enhancement roughly tracks that of cold dark matter is a robust 
result for small fi and large-scale loops. 

In the radial infall model for the Galaxy the magnitude of the cold dark matter enhance- 
ment is ~ 10^'^ at ~ 10 kpc. This quantitative degree of enhancement, although some- 
what model dependent, implies local loops with Gfi/c^ < 10~^^ are enhanced 10^ — 10^'^ 
with respect to the homogeneous universe. The degree of enhancement at the same 
position is also substantial for all /i less than the critical value for clustering; it exceeds 
103-6 for Gfi/c^ < 10-^0. 

The enhancement of the energy density of loops as a function of galactocentric distance 
for a range of possible /i is summarized in figure These loops are long-lived but 
ultimately transient residents of the galaxy. Experiments sensitive to the local Galactic 
population of loops, especially microlensing, will enjoy increased detection rates com- 
pared to homogeneous estimates. 

Small-scale loops produced via cusps do not cluster strongly because they do not live 
long enough for the universe's expansion to damp their initial relativistic motions. The 
galactic enhancement of loop (energy) density is ~ 10~^ — 10~^ that of cold dark matter 
depending upon the string tension (G/i/c^ < 10~^°) and some uncertainties in the velocity 
dispersion of the newly formed ultra-relativistic loops. At 10 kpc the net enhancement 
of bound loops is <1 — 10. Experiments searching for evidence of small-scale loops 
need to be sensitive to the homogeneous distribution throughout the universe; the local 
enhancement is small. 



-2- 



Keywords: String theory and cosmic strings, Strings and branes phenomenology. 



D-branes, Cosmology of Theories beyond the SM, Classical Theories of Gravity . 



Contents 



|T|. Introduction § 

[L1| String Tension |^ 

|2|. Executive Summary ^ 

1^. Halo Formation Model 



p Aim 

3^ Self-similar Radial Infall |^ 

37^ Model Specifics |^ 

^. Equations of Motion in Inhomogeneous FRW Cosmology |14 

|5|. Characterization of Orbits |17 



5.1 Radial Geodesies 0/ 



5.2 Rocket Effect 



^ Critical Acceleration 



§. Network Evolution ^ 

0| Status of Large and Small Loops ^ 



6.2 Birth Rate Density 34 



67^ Fragmentation and Cusp-Mediated Loop Formation 



0. Loop Clustering ^ 
|7.1| Probability of Capture 37 
|73 Halo Profiles ^ 



7.3 Truncation by Rocket 



7.4 Resuhs 



Current Halo Profile ^ 

Measures of the Loop Distribution ^ 

Point: Large Loops from Fragmentation Model [44 

Counterpoint: Small Loops by Cusp Formation pi 



1^. Conclusions 

1^. Equations of Motion in Inhomogeneous FRW 
|B|. Relating Center of Mass and FRW Accelerations 



- 1 - 



|C|. Variants of Figure 9 
|D| . Monte-Carlo Dynamics 



1. Introduction 

Three recent developments motivate an examination of the clustering of superstring loops 
on galactic scales. First, braneworld cosmological scenarios that provide a framework 
for inflation and the big bang, necessary ingredients in modern cosmology, generically 
produce string-like defects. Second, string theory allows and recent investigations of 
warped throat compactiflcation suggest that superstring tensions can be much smaller 
than the GUT scale: G/i/c^ << 10~®. Third, recent high resolution simulations of 
string networks suggest that ~ 10% of the energy available in horizon-crossing strings is 
transformed into loops within a few orders of magnitude of the scale of the horizon. 

In the original GUT scenarios, a phase transition at the GUT energy scale created 
string-like defects with tension G/i/c^ ~ 10^^ whose dynamical motions generated den- 
sity perturbations ultimately responsible for large scale structure formation ^ |^, 
Intercommutation chopped long, horizon-crossing strings into loops moving at relativis- 
tic speeds. These loops, distributed in the universe in an approximately homogeneous 
fashion, evaporated by gravitational wave emission within a few Hubble times. 

By contrast, the three advances lead to a qualitatively new and different picture for 
the fate of the loops. Large loop size and small string tension implies that loops survive 
for many expansion times. As such, they slow down by cosmic drag and fall into existing 
matter potential wells [|, ^. 

This paper shows that low tension string loops cluster and form a halo about the 
Galaxy. The enhancement relative to the homogeneous loop distribution is substantial 
both in terms of loop numbers and loop energy (total length). 

Superstring loops are roughly analogous to stellar objects in two respects: they have 
flnite lifetimes as local luminous sources of gravitons much like nuclear burning stars emit 
photons over a flxed main sequence lifetime, and they are massive, compact and optically 
dark much like the remnants formed in post main sequence evolution. It is interesting 
to consider how loops, like main sequence and post main sequence stars, may reveal 
themselves either directly by their intrinsic emissions or indirectly by altering photon 
propagation of background sources. Designing and planning searches for superstring 
loops will require a good understanding of the distribution of loops throughout the 
universe and especially our own backyard. 



-2- 



1.1 String Tension 



Inflation is an essential ingredient in modern cosmology]^, 0. In superstring theory a 
speciflc realization is brane inflation and the simplest example of brane inflation involves 
the interaction of a D3-brane moving toward a Z)3-brane sitting at the bottom of a 
warped throat [|l( 



ni III III, 0, IT3, Ig, The collision and annihilation 

of the brane pair initiates the hot big bang. Cosmic superstrings (F- and D-strings and 
their bound states) are produced and stretched to enormous scales [p2| , 



, gl, |2|, gl, 123, 

, |31|. After the epoch of inflation these superstrings evolve by the processes of 



28, 29 



intercommutation and gravitational wave emission to yield a scaling network in which 
there exists a stable relative distribution of long, horizon-crossing strings and sub-horizon 
loops I 



g, H, g, 0, |3|, 13, m, |9|, 11, ig, ig, |5|. 

The key property of a cosmic string is the tension /i, or in dimensionless terms, the 
string's characteristic gravitational potential G/i/c^. Theoretical understanding of the 
characteristic tension likely to emerge in a physically realistic string theory solution is far 
from complete. Initial estimates suggested 10~^^ ^ Gyu/c^ ^ 10~^ but recent analyses 



of multi-brane, multi-throat scenarios have effectively removed the lower bound [24, 26, 

mm- 

Empirical upper bounds on Gfi/c^ have been derived from null results for exper- 
iments involving lensing[46, ^ 50, 51, 52, |5^, gravitational wave background 



and bursts [0, |5|, |5|, |6|, |3, H, |5|, H, ||, H , pulsar timing H, ^ ||, Pand 

,1^ 



74, 75 



m 



cosmic microwave background radiation 
These may be generally summarized as follows: (1) Searches for signatures of optical 
lensing in flelds of background galaxies imply Gfi/c^ ^ 3 x 10^''. The analysis relies 
on the deflcit angle geometry of a string in spacetime and the accurate estimation of 
survey selection effects. (2) Modeling of the CMB power spectrum yields Gfi/c^ < 10~^. 
The limit is based on well-understood properties of large-scale string networks although 
the precise quantitative results are sensitive to unknown details of the spectra of string 
bound states and the probability of string-string interactions. (3) Pulsar timing stability 
gives Gfi/c^ < 10~^. The limit assumes that loops of near-horizon scale are created by 
the string network. 

In short, cosmic superstrings must have tensions substantially less than the original 
GUT-inspired strings and there is no known theoretical impediment to the magnitude of 
G^/c^ being either comparable to or much lower than the current observational upper 
limits. 

The lifetime r of a loop of size / to emission of gravitational radiation is, on dimen- 
sional grounds, r ~ Ic/ {Gfi). Smaller tension yields larger r. In a cosmology with power 
law growth in the scale factor, the scaling solution chops long strings into loops of size 
I ~ act, i.e. proportional to the size of the horizon. Here, the dimensionless constant 
a typifles a characteristic loop size from a broad, possibly multi-peaked, distribution of 
string loop scales. A key parameter is the number of expansion times before the loop 
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evaporates Ht ~ a/{G^/c^) where H is the Hubble constant. When Ht is not big the 
loops evaporate before the universe has expanded significantly. Such is the case for the 
usual GUT-inspired structure formation scenarios. 

The notions that the tension might be very small, G/^/c^ << 10^'', and that the loop 
size distribution include some objects comparable to the scale of the horizon (10~^ < a < 
10"^) yields qualitatively new cosmological features. If Ht » 1 the resultant string 
network will contain many old loops. Loops born with relativistic speeds are significantly 
slowed by cosmic drag before evaporating. A necessary condition for clustering on a 
galactic scale is that the loops damp to speeds less than the typical speeds within the 
galaxy: v/c < {y/c)gai ~ 10^'^. This is only possible for large Ht. 

This paper presents a detailed calculation of the infall of loops in a schematic model 
for growth and formation of a galactic scale cold dark matter perturbation in the presence 
of a scaling string network. ^ The strings satisfy the Nambu Goto equations of motion 
(hereafter NG strings). The loops accumulate and form a large halo similar in many 
respects to the Galaxy's gravitationally dominant dark matter halo but with several 
unique features: old captured loops decay by emission of gravitational radiation and 
are eventually ejected by the recoil associated with the anisotropy of their gravitational 
wave emission ("rocket effect") while new loops are added continually near the Galaxy's 
turn-around radius. At each epoch the galaxy is dressed with a long-lived halo of string 
loops. 

§2 outlines and summarizes the various calculations and the results. Modeling details 
follow in subsequent sections: §3 describes the halo formation model, §4 sets up the 
equations of motion for string loops in inhomogeneous Friedmann-Robertson- Walker 
(FRW) cosmology, §5 characterizes the free motion of loops, accelerated motion and 
derives an approximation for the critical acceleration that unbinds a loop, §6 describes 
the model for the string network, §7 characterizes the clustering for string loops generated 
by the network at a given time, with given size, and possessing given tension, and §8 
calculates the loops currently distributed within the Galaxy in two different models for 
loop formation. 

Finally, §9 discusses some of the implications and outlines future work. 



2. Executive Summary 

The background cosmological model is Einstein-de Sitter with a critical density of non- 
relativistic matter and scale factor a{t) oc t^^^. The growing gravitational perturbation 
is much smaller in size than the horizon and treated non-relativistically. Consider the 
spherically symmetric infall of cold, collisionless matter caused by introducing a small 



^ There have been numerous investigations of loop dynamics in homogeneous background! 54, 80, ^ 
|8^ , This work differs by fully accounting for the presence of growing gravitational perturbations 
which are ultimately responsible for the clustering. 
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overdense top hat perturbation. The turn-around radius is the point where the Hubble 
expansion exactly balances the infall velocity and sets the scale for the problem. Ma- 
terial well within the turn-around radius has had sufficient time to collapse, re-expand, 
recoUapse, and so forth. These motions ("bounces") rearrange the mass, kinetic energy 
and potential energy so as to produce a virialized structure on scales somewhat less than 
the turn-around radius. Material outside the turn-around radius has not yet had time 
to bounce. If much more material has bounced than was initially present in the top hat 
perturbation then the solution is self-similar, i.e. all the physical properties can be scaled 
from one time to another. The essential character of the self-similar infall model (i.e. the 
fully self-consistent density, velocity and gravitational potential of the collisionless cold 
dark matter which accretes and forms the non-linear bound object) has been spelled out 



[M ; |85| , |86| , |87[| . In this paper, the Galaxy's growing gravitational potential is described 
in terms of the gravitational potential ip of the self-similar infall model. 

The background metric for fiat FRW is ds"^ = —c^dt^ + a{tY{dx'^ + dy"^ + dz^). 
The perturbed metric is assumed to include only the scalar terms of the growing, non- 
relativistic gravitational potential ds^ = — {1 + 21^)0^ dt"^ + a{ty {1 + 2^) {dx'^ + dy'^ + dz^). 
Anisotropics in the stress tensor are ignored; this implies $ = There is no back- 
reaction of the strings on the spacetime. In the absence of gravitational wave emission, 
each string loop center of mass follows a geodesic in the spacetime; tidal effects on the 
loops are ignored. 

However, the loop does emit radiation and its impulse alters the loop trajectories. 
Intrinsic variations in the direction of momentum radiated within the loop center of mass 
frame would tend to average out the net impulse given to the loop. An assumption that 
no torques act constitutes the "worst case" for binding of loops to the Galaxy because 
the direction of the rocket force is influenced only by special and general relativistic 
effects which are always active not by additional, intrinsic variations in the loop center 
of mass frame. In this paper no internal or external torques act so that the direction of 
the rocket impulse in the center of mass frame behaves like the spin of a particle. The 
impulse direction is Fermi transported along the spacetime trajectory of the loop. 

The rates of gravitational wave emission of energy and momentum by oscillat- 
ing loops have been previously calculated for a sparse sample of loop configurations 



B^, 0, In the center of mass frame the approximate, period-averaged 
consequences are (1) a constant rate of change of length (or energy), dl/dt = V e{Gii/ c^)c, 
and (2) a constant degree of anisotropy which induces an impulse = Tp{Gfi/c^){c^/l). 
Based on the numerical studies, typically F^; ~ 50 and Fp ~ 10 but these may vary by 
a factor ~ 2 for individual loop configurations. ^ 



^This estimate for Tp may be systematically too big Large loops have had their cusps excised 
and the anisotropy of radiation emitted by the kinks that remain is smaller than would be generated 
by cusps. Hence, the force of recoil on large loops may be smaller than inferred from the calculated 
examples. The value used is conservative in giving the "worst case" scenario for binding of loops. 
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In summary, a loop undergoes non-geodesic motion in spacetime because of the 
rocket effect with Fermi transport of the impulse direction. In the center of mass frame, 
the loop shrinks according to a simple, approximate description and suffers a recoil based 
on a fixed degree of anisotropy which determines the magnitude of the non-gravitational 
acceleration (hereafter "the impulse"). 

When loops are chopped off from the string network they are typically moving at 
relativistic velocities. Initial conditions are drawn from a homogeneous distribution in 
space, with a relativistic center-of-mass velocity and a range of sizes. A choice of string 
tension /i and initial conditions for the loop (position, velocity, rocket direction and size 
at the time of formation ti) yields a trajectory in the background spacetime. 

The numerical results presented in this paper show how loops bind to the growing 
galactic perturbation. The most common scenario is that a loop born at an early time, 
slows down by cosmic drag, is overtaken by the turn-around radius, and accretes. The 
rocket effect is initially negligible. The resultant orbit is very radial passing back and 
forth through the galactic center with roughly fixed physical semi-major axis. Such a loop 
is easily identified as "bound to the perturbation." The physical scale of the captured 
orbit is fixed even as the perturbation continues to grow in size. Loops captured at early 
times end up near the center of the structure, ones added later at the periphery. 

As the loop shrinks, eventually, the acceleration of the rocket unbinds the loop 
from the perturbation. Detachment is rather sudden because the periodic motion in 
the potential averages the effect of the impulse, i.e. the orbit is adiabatically invariant. 
Escape follows when the rocket's force is large enough to alter the orbital parameters 
within a single orbit. 

The general trends can be understood by reference to the behavior of cold dark 
matter. Cold dark matter falls into the growing perturbation and creates a well-defined, 
universal profile with scale set by the turn-around radius. The cumulative probability 
distribution for a cold dark matter particle to be bound to the perturbation is a function 
of a single dimensionless parameter, the ratio of radius to turn-around radius. The 
analogous distribution for bound loops is more complicated since it depends on the 
epoch of formation, the size of the loop, and the string tension. 

Loops - formed at early epochs but with small enough /i that they have not yet 
neared the end of their lives - behave just like cold dark matter. They assume the same 
universal profile near the center while at large galactic radii the profile is truncated due 

to the rocket effect ("outer cutoff"). 

Dynamical complications ensue for loops formed at much earlier or much later times. 
For old loops that have begun to shrink substantially the importance of the rocket 
effect increases; consequently, the outer cutoff in the Galactic profile of loops becomes 
more important, i.e. it moves inward. Every loop is unbound before it completely 
evaporates because in that limit the acceleration from the rocket grows large. New 
loops, on the other hand, which are typically born with relativistic speeds have not yet 
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damped sufficiently to allow capture by the potential. A set of graphs illustrate the 
constraints on formation size, formation time and string tension. 

An integration over the rate of creation of loops implied by recent network simula- 
tions of Nambu Goto strings yields an estimate of today's loop proffie about the Galactic 
center. The main part of the bound population was created by fragmentation of horizon- 
scale strings into large sub-horizon loops. The loops are abundant and overdense with 
respect to the universe's average number and energy density of strings of all sorts. 

The bottom line results for large loops are presented in figures (number density 
of loops within the galaxy relative to the homogeneous value) and (energy density 
of loops relative to homogeneous value). These figures show a substantial degree of 
enhancement in both measures that depends upon string tension and spans a large 
interval of galactocentric radii. Such proffies motivate more accurate calculations for 
gravitational wave, pulsar timing and microlensing experiments hunting for evidence of 
loops within the Galaxy. 

The radial distribution of loops within the Galaxy is weighted to the center. The 
size distribution of loops is weighted to small scales with a cutoff corresponding to a loop 
with evaporation lifetime equal to the age of the universe. 

By contrast, cusp-generated small loops fail to bind to the Galaxy. This is not 
surprising given their ultra-relativistic initial motions and their reduced lifetimes. The 
bottom line results for small loops are presented in figures ^ (number density) and 
(energy density). Little enhancement is observed within the range of radii at which the 
radial infall model is applicable. 

These results are broadly suggestive that clustering of large loops will play an im- 
portant role in setting microlensing rates and may also increase the effective sensitivity 
of gravitational wave and pulsar timing experiments. On the other hand, experiments 
sensitive to small loops will not benefit from significant local enhancement. 



3. Halo Formation Model 

3.1 Aim 

The role of the halo formation model in this paper is to provide a dynamical background 
for the motion and eventual capture of the string loops generated by the network. The 
density and potential of the cold dark matter are determined in a self-consistent fashion 
and the string loops moved as test particles in the potential. The distribution of cold 
dark matter particles and string loops can then be compared. 

3.2 Self-similar Radial Infall 



Refs. [Q, |8^, ||86|, and |^ have analyzed the spherically symmetric infall of cold. 



collisionless matter onto small-scale density enhancements in an Einstein-de Sitter uni- 
verse. The solution is self-similar i.e. the form and appearance at any time is fixed when 
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scaled to a characteristic physical length. The infall yields power law halo density dis- 
tribution p oc r~^'^^ and rotation curves oc r~^/^. The model is physically self-consistent 
and simple enough that many difficult aspects of the cosmology plus network evolution 
can be handled precisely. 

The model solution is exact given the assumptions but the model must be regarded as 
a schematic description of the Galaxy. Its basic shortcomings in comparison to a realistic 
treatment of ACDM cosmology are the following: (1) it ignores the initial spectrum of 
perturbations which span a great range of length scales and it suppresses the generic 
asymmetry of a typical perturbation that grows to encompass a galaxy scale mass, (2) 
it is only applicable at times after equipartition and before late-time acceleration. 

Of these, the more significant issue is the ffist. Initial conditions drawn from a 
CDM spectrum generate a hierarchy of mergers not a monolithic infall [p^ , |95| , p7| . 



Today objects like our Galaxy have dark matter density profiles that are non-power 
law The density p varies like 1/r at small radii (successive mergers of small dense 
objects) and at large radii (truncation of infall). The radial infall model does not 
capture the behavior at either extreme: it should be adequate on scales on which the 
rotation curve is observed to be fiat, 3 < r < 30 kpc, and it may be adequate out to 
distances where the curve is traditionally assumed or inferred to be fiat, e.g. ~ 100 
kpc[^]. The motivation for its use here is that a comparison of the loop and cold dark 
matter distributions calculated in the same, self-consistent time-dependent potential, 
should yield valid conclusions of greater generality than might be suggested by a strict 
comparison of actual to modeled Galactic proffies. 

Of course, a oc t^^^ prior to the time of equipartition t^q and the perturbation begins 
to grow only for t > teq- It is incorrect to use the self-similar form for the potential at 
early times. Nonetheless, this paper employs it focusing on clustering after equipartition. 
Since the dynamics of loops formed before equipartition are an essential part of the story 
to be set forth one might worry that this presents an additional important shortcom- 
ing. In the current model, the turn-around radius at t^q is ~ 15 pc; the radial profile 
of cold dark matter and of loops at such small scales would certainly be inaccurately 
represented even if the perturbation were exactly spherical. (Which it is not - as already 
indicated the smallest radius at which the radial infall model applies is much larger.) On 
larger scales the distribution of cold dark matter and loops is not adversely impacted 
since the potential is basically flat while cosmic drag and the rocket effect both operate 
independently of the potential. A loop slows down first and then binds to the growing 
perturbation when the turn- around radius reaches it. This typically occurs at t ^ t^q. At 
times t > teq, the cold dark matter density and gravitational potential quickly asymptote 
to the self-similar form. Computed properties today are uninfiuenced.^ 

Finally, the late-time acceleration of the universe alters the behavior of the turn- 



^An independent issue relates to the string network evolution and the distinction between scaling 
solutions before and after equipartition. This is discussed later. 



- 8 - 



around radius. This certainly changes how loops might be added to the outermost 
periphery of the Galaxy. Since M31 will turn out to be closer than turn-around radius 
formally inferred today, it is clear that quantitative agreement at such large distances is 
not expected. 

Future work will treat structure growth in realistic simulations of the ACDM paradigm. 
3.3 Model Specifics 

This exposition follows the mathematical description given by Ref. [^. At time ti 
consider a small top hat density enhancement 6i = {Sp/pH)i which extends from the 
origin to physical radius Ri; assume the Hubble constant Hi is independent of radius. 
All the material in the universe is bound by the presence of the excess material and is 
destined to fall towards it. 

The motion of shells depends upon the initial position rj, the initial velocity rjifj 
and the mass interior to the shell. Write the initial mean density within a radius r,- in 

shell, it satisfies the parametric equations 



2 (a - 1) 



2H, (a - 1) 



1-COS0) (3.1) 
-sine). (3.2) 



3/2 



The initial state at t = has parameter ^ = 0. To lowest order in << 1 it turns 
around {6 = tt) at 

r. = ^ (3.3) 

1/3 

Eliminating the occurrences of in favor of 6^ Ri yields the turn-around radius at any 
given time t as 

n. = 6y'R.(^^y\ (3.5) 

The fact that Vfa involves the product of a single combination of variables and the 
characteristic power law t^^^ is a significant simplification. The turn-around radius is a 
suitable length scale with which to non-dimensionalize the problem. The dimensionless 
length scale A is a pure function of the parameter 6: 

, rjt) sm^(l) .36. 
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A shell starting at a large distance from the top hat has large A. As the turn-around 
radius increases and the shell falls under the influence of the perturbation A decreases. 

The initial motion is orderly. Before shell-crossing occurs the dimensionless velocity 
V, dimensionless mass and dimensionless density p are explicit functions of 6: 

i; ^ ^^^^ ^^^^ (Q ■ ml/9 

_ 3M _9A Mg-singf 

AnpnityUtr 16 sin^ (|) ^^'^^ 



p _ 1 (K^-sin^))' 



PHit) 1 + 3x sin^ (f ) 



(3.9) 



3V 

X- = 1 - ^ (3.10) 

and pnit) is the density of the background model. 

The initial state corresponds to ^ ^ 0; the asymptotic forms in this limit are 



A - 


4^2/3 


(3.11) 




2A 


(3.12) 


V - 


Y 


M - 


A^ 


(3.13) 


P - 


1. 


(3.14) 



The excess mass associated with the perturbation 

is used to infer the potential that will make its appearance in the metric in the next 
section. 

Eventually, the infalling shell meets shells that have already bounced. The first 
crossing for a shell occurs at Ai = 0.364. The dimensionless expressions above are exact 
for A > Ai. Once crossings begin for the shell of interest, the mass interior to it varies 
and must be calculated to find the acceleration and trajectory. Ref. has solved the 
problem numerically and tabulated A^(A) for 0.02 < A < Ai and also provided a power 
law approximation for A << 1. This paper uses a combination of asymptotic, tabulated 
and analytic expressions to describe A^(A) and Ai^W over the complete range of A.'' 
The dimensionless cumulative mass is shown in figure |l|. 



■^The forms for A << 1 follow the scaling given in Ref. |Q with coefhcients adjusted to fit the last 
tabulated values; the numerical coefficients differs at the 5% level from limiting analytic expressions 
given elsewhere in the same paper. 
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The dimensionless radial component for the force and potential are deduced from 



f 



fGpit)rtait) 



A2 



'-d\. 



(3.16) 
(3.17) 



The fact that the force and potential, which are intrinsically functions of four spacetime 
variables, are compactly represented in terms of the one-dimensional functions of the 
dimensionless radius A = r/rta{t) is a great simplification. Test particle motion in 
the vicinity of the perturbation depends upon the these functions. The dimensionless 
potential is shown in figure 




Figure 1: The cumulative mass bound to the perturbation as a function of distance from the 
center; A = r/rta is the scaled radial coordinate; the mass is expressed in units of A'irp{t)rf^/3 
where p{t) is the background (critical) density. The curve is continuous with changes in slope 
(density discontinuities) at shell-crossings. The upper curve is the entire mass including that 
contributed by the mean background; the lower curve is the excess with respect to the back- 
ground. 



To apply the infall model to the Galaxy today the turn-around radius must be 
specified. Assume that the well-virialized part of the halo has a physical scale today 
Rg = 150 kpc. The current age is fixed at the concordance value to = 1.37 x 10^° yrs 
Let the epoch for the turn-around of the material at Rg be ti and let the turn- 



100 



around radius at that time be rta,i- Figure || illustrates the history of a particle as it 
passes back and forth through the center. Its apocenter approaches an asymptotic value 
of ~ 0.8rta,i. Between 4 and 5 passages, the apocenter has shrunk from rta,i at time ti 
to 0.845rta 1 at time 8ti. With some arbitrariness, identify Rg = 0.845rta,i and to ~ 8^i- 
Since r^a oc t^^^ this implies the turn-around radius today is rtao = 1-1 Mpc. 
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Figure 2: The gravitational potential associated with the excess mass as a function of distance 
from the center; tf^x is the gravitational potential measured in units of 47rG/9(t)r^„/3 (with value 
zero at spatial infinity) and A = r/rta- 



All the quantitative properties of the self-similar description now follow. Assuming 
that radial infall continues to the present, at the Sun's galactocentric radius today (8.5 
kpc) the rotation velocity is 225 km s~^ and the interior mass is IO^^Mq. The model 
has total mass 7.5 x lO^^M© within Rg today. Turnaround at ti/to ~ 0.12 occurs before 
the switch from power law to exponential expansion in the concordance ACDM model 
(^a/^o ~ 0.73) so the late-time deviations from Einstein-de Sitter should be relatively 
unimportant. If infall were to cease at t > ti the net change in mass within the solar 
circle is estimated to be ~ 4% and within Rg to be ~ 40%. 

This completes the specification of the time-dependent Galactic model. 
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Figure 3: The motion of a single particle in the collisionless infall solution. Time is in units 
of the initial turn-around time ti, physical radius in units of the initial turn-around radius rta,i- 
The curve is a fit to numerical results (Table 6 in |]87|| ). 
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4. Equations of Motion in Inhomogeneous FRW Cosmology 



The unperturbed background model is flat Friedmann- Robertson- Walker (FRW). Con- 
sider a frame in which the universe appears isotropic and let = (t, x*) where t is the 
global time coordinate and are the global comoving spatial coordinates. The scale 
factor is a{t). Henceforth, adopt units with c = 1. 
The FRW metric with scalar perturbations is 



(-{1 + 2^) \ 
a2(l + 2$) 

a2(l + 2$) 

V a2(l + 2<l>)/ 



(4.1) 



where il){x'^) and $(x") represent the effect of inhomogeneities. Bound structures have 
^ < 0. 

The perturbation potentials are very non-relativistic (rotation velocity ~ 225 km/s 
for our Galaxy implies v ~ 10~^). The Galaxy is assumed to be at rest in the preferred 
FRW frame. Only small errors are made by equating the local time and space coordinates 
introduced in the previous section to the global FRW coordinates {t,ax^). Assume also 
$ = —tp which is suitable for small quadrupolar components to the stress energy tensor 
for the cold dark matter particles. 

In the FRW frame, let a loop's velocity be V^. If there were no rocket effect the loop 
would follow a geodesic through spacetime (ignoring tidal effects). However, the loop 
does emit radiation and its impulse alters the loop trajectory. Let the rocket's 4-impulse 
be arN"" where is the time-varying magnitude of the impulse. 

Consider, first, the direction of the rocket's impulse. In the center of mass frame, 
(1) loops formed in cosmological fragmentation scenarios generally possess net angular 
momentum ||101|| , (2) the angular momentum radiated over an oscillation period lies par- 



allel to the angular momentum of the loop[^, and (3) the momentum radiated over an 
oscillation period lies in a direction generally different than that of the angular momen- 
tum of the loop. Item (2) implies that loops spin down in a relatively simple manner, 
however, item (3) suggests that the net gravitational force does not act on the center of 
mass of the loop. The emission of angular momentum has been studied only for loops 
of the simplest complexity. Though there is no evidence to date, more complex loops 
might experience more complex dynamics. 



Ref. ||80| treats the loop as a relativistic gyroscope and concludes from a dimensional 
argument that the timescale for a single precession cycle is ~ r, i.e. comparable to the 
loop lifetime. Since the momentum impulse is not along the angular momentum direction 
this argument suggests but does not prove that the rocket direction is fixed for the life 
of the loop. On the other hand, if one treats the loop as a solid body |^ the precession 



time is considerably shorter ~ y^Gfi/c^r. A suitable gravitational wave back-reaction 
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calculation that would definitively address how the direction of momentum impulse varies 
is unavailable ( ||102|| did not investigate precession and ||103|| studied a symmetric loop 
that did not radiate momentum). 

Intrinsic variations of the rocket direction in the loop center of mass frame would 
tend to average out the net impulse given to the loop. A fixed direction gives the "most 
effective" rocket and presents the "worst case" for binding of loops to the Galaxy. 

Assume (1) the impulse in the loop center of mass frame lies along a fixed direction 
and (2) no torques act in the loop center of mass frame. Then A^" is simply the Fermi- 
transported impulse direction of the loop. The normalizations are V^Va = —1 and 
N"Na = 1 and orthogonality is N^Va = 0. The equations of motion are 

-j^ + r^^y^V^^ = a.AT" (4.2) 

-j^ + r^^AT^T/^ = a,\/" (4.3) 

for proper time r. For the numerical solution in the FRW frame, the 4- vectors for velocity 
and for the internally generated impulse direction are parameterized 



= I 1 (4.4) 



m = ( ±^^, ""^ ) (4.5) 

where v"^ = QijV^V^ , = gijN^N^, and -0* and n* are 3D-orthonormal unit vectors 
[v ■ V = Y^-v'^v'^ = 1 and n ■ n = 1).^ The equations for dx^/dt, dv/dt, dv/dt, dn/dt 
and dh/dt are expressed using the global FRW time as the independent coordinate. 
These equations are applicable to loops with the whole range of possible velocities from 
extremely relativistic to non-relativistic. The explicit form is given in the Appendix ^ 
Let the total loop energy be E in the FRW frame and, following custom, denote E/fi 
as length /. For clarity, explicitly label quantities in the string's center of mass frame 

with "z". The infinitesimal length (i.e. energy) is dl(z) = da\J (dz/da)"^ /(I — ^) where 
z = z{a, z^) is the parametric expression for the string; / = V^^z) and dt = V^dz^. 

In the loop center of mass frame the rate of energy loss and the magnitude of the 
impulse are very simple 



dl(z) 



-TeGix (4.6) 

rp^. (4.7) 



^With this parametrization, an FRW observer sees an energy per mass vl + i? and a momentum 
per mass v. In terms of relativistic kinematic variables v = jP. In this paper, v is called "velocity" 
when the regime is non-relativistic and "momentum-per-mass" for more generality. 
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The loop lifetime is a fixed increment of time in the center of mass frame. The 4-impulse 
in the center of mass frame is a"^-) = (0, a*^^-)) = ar(0,n*^^-|) where n*^^^ is a unit vector in 
the direction of the impulse and is the magnitude of the impulse. Since is a scalar, 
a'^act = a,'^^-^a,{z)a = ci^-, write the 4-impulse a^-A^" in the FRW frame. Since the initial 
loop configuration z{a, 2;°) determines n*^^^ and no torques operate (by assumption) it 
is most convenient to find the initial N"" in FRW frame (Appendix ^ and use Fermi 
transport to determine its subsequent evolution in that frame. 
In the FRW frame 

dl I dV^ ^ ^ N 
- TeG^i (4.8) 



dt 

d f l\ Te 



dt\ar TpV^' 



(4.9) 



After a^A^" is initially set the entire calculation can be carried out in the FRW frame 
using Fermi transport for A^"" and the above equation for a^. 

There are a variety of non-trivial frame transformation effects that operate in this 
schematic description of loop evolution. Ignoring the momentum impulse of the rocket 
and the inhomogeneous potential, in the FRW frame a loop born with length Zj with 
center of mass motion = y^l + vf at time has length and center of mass momentum 

TEGfi r ; r/^' x^l^dx , 




(4.11) 

at scale factor a = ai{t/tiY . The initial loop size that just evaporates at scale a/a^ is 
explicitly given by setting the expression within the parenthesis to 0. It is clear that 
complete evaporation occurs in a finite FRW time. 

The time a loop lives is slightly different than the above result in a homogeneous 
universe because of the ever increasing importance of the rocket effect. The loop length 
still vanishes in a finite time. Let the time until evaporation be At = tuje — t and 
the ratio of momentum-to-energy loss (3 = Tp/Te- Then the length, acceleration and 
momentum parameter vary asymptotically 

/ oc At (4.12) 

a, oc At-i/(i-^) (4.13) 
V oc At-^/(i-^) (4.14) 

The difference between the approximate and exact lifetimes is only ~ 2.3% (F^; = 50, 
Tp = 10, Vi = 0.1) which will be ignored in subsequent discussion. The acceleration and 
the momentum-per-mass both diverge as the evaporation proceeds to completion. 
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5. Characterization of Orbits 
5.1 Radial Geodesies 

The first calculations illustrate some basic kinematic features for objects whose initial 
velocity is very different from Hubble flow. Consider purely radial geodesies and ignore 
the rocket effect. Fix the magnitude of the initial radial velocity to be fj = 0.1 (typical 
of the largest loops chopped off from the horizon crossing strings) at the initial time 
ti/to = 10~^. The results that follow are for the full, relativistic equations of motion. 
One calculation differs from another only in terms of the initial position of the loop with 
respect to the spherical center. 

Figure | shows the comoving trajectories of a set of loops as a function of logt/to- 
The different colors label different initial positions and ingoing and outgoing velocities. 
The velocity v is damped by the many decades of expansion (in the absence of a varying 
potential v{t) ~ Viai/a{t)), a fact made qualitatively clear from the flattening of all the 
curves at early times. Once a loop's motion has been damped, it behaves for all practical 
purposes like a cold dark matter particle at the position to which it has moved. The 
comoving coordinate is nearly but not exactly static because every zero-velocity object 
is bound to the excess central mass of the perturbation. In the absence of the rocket 
effect, each loop eventually turns around, the comoving coordinate retreats and the loop 
oscillates back and forth through the perturbation center. To avoid clutter, only the first 
few bounces of each loop are plotted. The color of the line allows tracing the epoch of 
turn-around and recollapse for a loop to a given initial position. The black line is the 
comoving turn-around radius in the radial infall model. 

Figures |^ and ^ are blow ups in comoving and physical coordinates respectively 
for 10"'^ < t/to < 1. They show the first few bounces after capture of the loops by 
the perturbation. These figures as well as the previous one demonstrate that early (late) 
turn-around implies small (large) semi-major axes just as is true for the cold dark matter 
particles in the radial infall model. Figure ^ illustrates that loops starting in different 
regions of space with different initial velocities (red and blue lines) can end up with 
nearly identical accretion orbits. This is simply the shuffling in position that occurs 
during the time it takes for cosmic drag to operate. As a side note, the absence of red 
lines is a consequence of the limited range of initial radii sampled; had larger offsets 
been plotted such lines would be present throughout the figure. Figure |^ makes it clear 
that there exist loops that "turn around" at the same space time locations as cold dark 
matter particles do. 

If a loop is young then cosmic expansion may not have had sufficient opportunity to 
damp its velocity to allow accretion onto the growing perturbation at a physical radius 
of interest. A simple estimate of how small ti/to must be for a loop to capture at time 
t is given as follows. Cosmic drag implies the initial velocity decays like v{t) = Viai/a{t) 
(flat potentials). A necessary condition for capture is that v{t) must be less than the 
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Figure 4: Comoving position as a function of logt/to- Radial geodesies (no rocket effect) 
have been integrated from a set of initial positions with fixed initial radial velocity Vi = 0.1 
(either inward or outward) at time ti/tQ = 10~^. Distance is measured in units of today's 
comoving turn-around distance. The black line shows the comoving turn-around radius. 



escape velocity from the potential. However, a generally more restrictive condition is that 
capture requires v{t) < H{t)rta{t). Since the turn-around radius rtait) = rtait, 







the initial time is constrained to be tj/to < ( " ] (^] • For example, a loop 

with initial velocity Vi = 0.1 can be captured today if tj/to ~ 7.5 x 10~^ and will have 

a physical orbit ~ rta{to)- Physical radius and time of capture are inherently linked in 

the similarity solution. A loop with physical orbit r < rta{to) must be accreted at earlier 

time t/to = (r/rja(to))^^^; the initial time of formation of that loop is constrained to be 

fH r it A^/^ f ■ \ ^^^^^ 
ti/io < ( ° ] irTlfy)) ■ "^"^^ example, for r = 30 kpc, sufficient cosmic drag 

requires the loop be formed at ti/to < 2.5 x 10~^ and captured at t/to ~ 1.7 x 10"^. 

Figure |^ shows the comoving trajectories of a set of loops born at ti/to = 10~^ (the 

ordinate is greatly expanded compared to previous figures). The slope for orbits far away 

from the perturbation indicates that cosmic drag has not yet brought the loops to rest. 

This impedes capture. Figure ^, a detailed view near the origin, shows that when it does 

occur it does so at large physical separation. 
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Figure 5: Detailed view of the trajectories in figure ^. The different colored lines indicate 
different initial conditions. Only the first few bounces through the center are plotted. Con- 
sequently, the increase in amplitude with time highlights a sequence of loops falling into the 
perturbation having systematically larger orbits not that an individual loop's bounces grow in 
amplitude. 
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Figure 6: Detailed view of physical separation for trajectories in figure ^. Only the first few 
bounces through the center are plotted. In physical coordinates, the bounces occur with nearly 
fixed amplitude. Same comments as for figure |5| apply. 
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Figure 8: Same as figure |^ but in physical coordinates. 
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5.2 Rocket Effect 

The rocket effect can inhibit binding of a loop to the growing perturbation and can 
unbind a previously captured loop. This section begins with analytic estimates and 
follows with full numerical calculations. 

Begin by considering homogeneous FRW with rocket direction aligned or anti-aligned 
with initial velocity. The equations of motion (Appendix 0) reduce to 

dv ^ nar a . , 

— = 0—. — V- (5.1) 

dn var nv'^a . , 

dar _ TeciI , . 

"dt ~ Fp^IT^ ^ ■ ' 

6 = n-v = ±l (5.4) 

with initial velocity Vi and loop length U = a/ Hi = (3/2)ati in the matter-dominated, 
FRW frame at time ti. The initial rocket impulse is a^^i = TpG^^Jl + and the FRW 
direction vector is ni = \/l + vf. The equation for n may be solved and its occurrences 
eliminated. For non-relativistic motions, the equation for a-r may also be integrated 
explicitly to give 

— = 0— V- (5.5) 

dt T E Ti + ti — t a 

n = ^ (5.6) 

Here is the characteristic time for the loop to evaporate by gravitational wave radiation. 
li Ti + ti » t then the first term is approximately constant and f is a sum of terms 
proportional to t and The limit of interest is Hr » 1. Specifically, if 

K-^ = {5TE/i2Tp))viHiTi = {5aVi/{2TpGfi)) » 1 then 

2/3 

vc^Vi{{"-^] +6'^] (5.7) 





Capture at radius r of the accelerated trajectory for a loop with aligned velocity and 
rocket impulse (6 = 1) requires 

Horta.o f tQ\^^^ f r \ _ ^ 
'^^ . N 15/8 / N 5/3 y^-^' 



^rta,0 J V / 

by the same line of argument given in the previous section. The numerator must be 
positive and implies an upper limit on the time of formation 

\ Vi J \rta,oJ to 
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identical to that in the previous section. This is the effect of cosmic drag and does not 
depend upon /i. For loops formed at earlier times, capture requires 



(I^)-^ (-0) 

a result that highlights the importance of the rocket effect. 

The condition that the loop not evaporate by the current epoch requires 

G^<^!i. (5.11) 
2Te to 

For a given epoch of formation, the upper limit on tension is rocket-related for r / rj^ q > 
{5T E / {'iT p))^/^ (Hortafl)^^^ ~ 7.6 x 10~^ or r > 8.5 kpc and age-related at smaller radii. 

The maximum tension that permits capture (the intersection of limits implied by 
equations and ^.10| or this is the rightmost section of the triangular region 



formed by green and turquoise lines in figure ^ ) is 

2Fp \ Vi J \ r J y 5Te \rta,oJ J 

Above this critical Gyu/c^, loops do not cluster at scale r within the Galaxy. Curiously, 
the critical G/i has a maximum close to the Sun's galactocentric position though the 
variation from r = 3 — 100 kpc is only about 2.7 (and with all other parameters fixed). 

All captured loops are eventually stripped from the galaxy by the rocket effect. A 
detailed discussion of how removal proceeds is given in the next section. The result is 
that ttr > 0.3|VV'| leads to detachment. For fixed loop size, a loop is retained until the 
current epoch if 

0.15a M^U 
Gj^ < — — Hovtafi—^— (5.14) 
Tp to 

For small A the asymptotic form for the mass distribution is Ai^ ~ 10. 5A^/^, so that 
Aix/ oc A"^/^ just as in eq. |5.10| . This retention criterion is very similar to the capture 



criterion but quantitatively a bit stricter (the constant for retention is ~ 0.6 that of the 
capture). However, the physical interpretations are very different. Capture is a statement 
about the properties of the loop and turn-around radius at early times whereas retention 
concerns all later times up to the current epoch. Because of the self-similar evolution 
of the perturbation both criteria vary with to identically for fixed loop size. As the loop 
shrinks, the retention criterion becomes more and more strict. In summary, the retention 
criterion is only a bit stricter for fixed loop size but becomes far stricter once the loop 
size begins to change. 
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Figure |] provides a graphical summary of some of the analytic results and a point 
of reference throughout this section which describes additional calculations.^ 

The above analysis is now supplemented by a study of the rocket's effect on loop 
dynamics via a sequence of trajectory calculations of increasing complexity. The nu- 
merical investigation proceeds along the following lines. First, locate the specific radial 
trajectories that give rise to orbits of fixed physical size today for /i = 0, i.e. without the 
accelerative force. Second, starting from the inferred initial conditions, repeat the cal- 
culation of the trajectories including the effect of gravitational wave recoil. The rocket 
effect is slight at first but dominates by the end of the loop lifetime. Comparison of 
different fi yields a criterion for retaining a loop at the radius of interest today. Figure 
H includes some of the numerical results. 

Begin by choosing a time for the birth of the loop tj/to (10^^ to 1) and fixing the 
initial radial velocity {vi = 0.1). Then find the initial position that ends with an orbit 
of the desired apocenter at the current epoch (two specific cases, 10 and 30 kpc, are 
considered; the apocenter is estimated based on the last 2 extrema of the orbit prior 
to to). This is a boundary value problem for the equations of motion and is solved by 
numerical iteration. 

For each case, two qualitatively distinct orbital solutions for small tj were found; 
no attempt was made to find all solutions. As ti increased the initial conditions for 
the individual solutions converged and for t greater than a critical value no suitable 
initial conditions could be found. This result is consistent with the order-of-magnitude 
argument given above. The variation of the initial position with tj/to of the two branches 



is systematic and shown in figure [10 



Both branches of 30 kpc orbits are over plotted in comoving coordinates in |TT]; they 
begin at different times and different velocities but note that all converge to similar os- 
cillatory solutions. These solutions are the "baseline" solutions which are now perturbed 
by the rocket effect. 

Starting from the initial conditions inferred above, repeat the calculations with non- 
zero fi and a set of random orientations for the rocket in the loop center of mass frame. 
Varying fi delimits the transition from bound to unbound orbits at the current epoch. 
Small (large) n implies weak (strong) acceleration. The transition refers to a specific 
time, t = to, as it is clear that eventually, for long enough integrations, the orbits of all 
evaporating loops are unbound. 

The initial position and velocity in figures p!^ - |l^ are all identical, i.e. that of a loop 
born at tj/to = 10^^ which formed a bound 30 kpc orbits in the absence of the rocket 
effect. A sequence of calculations with increasing log^g/^ — —12.8, —12.7 and —12.6 and 
random rocket orientations unbinds an increasing fraction of orbits. Each figure includes 
the original unperturbed orbit for comparison. 

Besides the fact that the transition occurs over a fairly narrow range in /i, the figures 



thank Xavier Siemens for his version of this figure. 
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Figure 9: Bounds on formation time and string tension for a loop with initial velocity 
Vi = 0.1 to be captured at physical radius 30 kpc. (1) Upper bounds on the formation time 
ti/to are given by the horizontal lines. The condition that cosmic drag lower the velocity to 
less than the circular rotation velocity today is given by the red line. The more stringent 
condition that capture occur at 30 kpc is given by the turquoise line. (2) Upper bounds on 
the string tension Gfi/c^ are given by the diagonal lines. The initial loop size is a/ Hi and 
a = 0.1. The condition that the loop be younger than its gravitational wave decay timescale 
is given by the red line. The more stringent condition that the loop not be accelerated out of 
the Galaxy is given by the green line (x > Xcrit)- The geometric symbols summarize results 
from numerical experiments examining the outcome today {t = to) for groups of 10-20 loops 
captured at 30 kpc with slightly different string tensions: stars = all loops bound, boxes = 
all loops ejected, triangles = some bound and some ejected (for clarity the points are slightly 
offset in the vertical but not the horizontal direction). (3) The triangular region encompasses 
string tensions and formation times giving bound loops at 30 kpc for Vi = 0.1 and a = 0.1 in 
a radial infall model of the Galaxy. The critical value of fi below which clustering is possible 
is determined by the upper right hand corner of the green and turquoise lines. Lowering Vi 
raises the limit on ti/to (horizontal lines moves upward); lowering a shifts the bound to smaller 
fi (diagonal lines moves leftward). Shifting the loop orbital scale to smaller values (say the 
solar position) requires earlier formation times (horizontal lines shifts down) and allows larger 
fj, (green line shifts to the right but is limited by the red line which is fixed). 

also illustrate that the apocenter does not significantly change until just before the orbit 
is actually destroyed. 

A summary of the results for a range of /i and ti/to for bound 30 kpc orbits appears 
in figure ^ The geometric symbols characterize the outcomes for groups of orbits: stars, 
boxes and triangles mean "all bound", "all ejected" and "a mixture" of both types, 
respectively. The green diagonal line is an analytic estimate for the critical n based on 
non-relativistic calculations in the following section which imply that a loop is detached 
when the magnitude of its rocket impulse satisfies ^ 0.3|V'?/'|. The transition from 
bound to unbound is abrupt (a factor of 2 in /i encompasses its entire scope) stemming 
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Figure 10: Initial comoving coordinate positions (ordinate) at a set of different formation 
epochs (abscissa) for loops that evolve to give orbits of fixed physical size today. The green 
(blue) points are the initial conditions that give 30 (10) kpc orbits. The two branches are two 
qualitatively different solutions to the boundary value problem described in the text. The loop 
in the lower branch passes through the perturbation center, slows down and is overtaken by 
the turn-around radius. The one in the upper branch begins moving outward and is overtaken 
by the turn- around radius. 
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Figure 11: The different lines illustrate the comoving radial coordinate of loops formed 
at different times and yielding 30 kpc orbits today. These are the explicit solutions to the 
boundary-value problem whose initial conditions are displayed in Figure |l^; they correspond 
to the green upper and lower branches. 



from the adiabatic invariance associated with averaging weak forces over periodic orbits. 

This figure summarizes the main physical constraints for binding and residency of 
loops in the Galactic halo. There are upper bounds on the formation time tj/to given 
by the horizontal lines. The condition that cosmic drag lower the velocity to less than 
the circular rotation velocity is given by the red line. The more stringent condition that 
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Figure 12: Initial time U/to = 10^^, initial velocity Vi = 0.1, fi = 10~"^^'^. Eight trajectories 
with randomly chosen initial momentum directions and the original unperturbed trajectory 
(golden yellow color) are plotted 
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Figure 13: Same as ^ but /Lt = lO"^^'^. One trajectory is detached. 

capture occur is given by the turquoise line. There are upper bounds on the string tension 
given by the diagonal lines. The condition that the loop be younger than its gravitational 
wave decay timescale is given by the red line. The more stringent condition that the 
loop not be accelerated out of the Galaxy is given by the green hne. The critical Gfi/c^ 
below which capture and retention is possible is given by the intersection of the green 
and turquoise lines. Appendix y includes a graph which replaces the straight lines with 
numerically determined conditions for aligned and anti-aligned rocket orientations. 

The triangular region encompasses string tensions and formation times for loops that 
are bound to our Galaxy today. The specifics of this figure refer to loops at radius 30 
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Figure 14: Same as but /i = 10~^^'^. All the perturbed trajectories escape. About half 
the trajectories are captured and later detached while the rest avoid capture completely. The 
unperturbed trajectory is golden yellow. 



kpc, with initial velocity Vi = 0.1, and initial length a/Hi where a = 0.1 in a radial infall 
model of the Galaxy. The manner in which the constraints vary is briefly discussed in 
the caption. 
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5.3 Critical Acceleration 

It is well known from classical mechanics that the action of a simple harmonic oscillator 
with frequency uj is an adiabatic invariant. Perturbative driving forces having intrinsic 
frequency Vt such that Vt « uj produce exponentially small changes in the action or 
energy. Here, perturbative means that the magnitude of the forcing is small, i.e. the 
instantaneous change to the coordinate is first order. The integrated change of a first 
order quantity over a full period is very small. 

The rocket acts on the loop's orbit within the Galaxy. For loops which have slowed 
enough to bind to the Galaxy, the Fermi transport of the impulse direction (an analog of 
Thomas spin precession) has frequency Vt ~ v'^uj ~ 10~^u;. Effectively, the acceleration 
is in a fixed direction with magnitude governed by the decrease in length of the loop. 
The condition for escape is equivalent to the breakdown in adiabatic invariance of the 
oscillator that occurs when the force grows sufficiently large to become non-perturbative. 

The main complication for an actual orbit is that the potential is not separable and 
several incommensurate tu's exist (radial and angular frequencies), so that one cannot 
solely focus on the motion in the coordinate direction defined by the impulse. However, it 
is straightforward to investigate a simple, non-relativistic model having all the essential 
features and to infer an approximate criterion for the transition from perturbative to 
non-perturbative motion. Consider an acceleration law of the form 

— * 

a = ""^ ^ +k. (5.15) 

(r^ _l_ j-T^ 2 

The interpretation here is that the first term is the acceleration due to the galactic 
matter distribution and the second term is the internal acceleration due to the rocket. 
The constants are core radius Tc, internal acceleration k and galactic acceleration power 
law n. For a Keplerian potential n = 2, for the radial infall model n = 5/4 and for 
galactic potentials the range of interest is 1 < n < 2. To non-dimensionalize, express 
lengths in units of Tc- 

For a numerical investigation, first choose the initial radial displacement of a zero- 
velocity particle and the size of the internal acceleration, and sample random choices of 
direction k. The unperturbed orbit would remain radial if not for the internal acceleration 
which drives it away from that limit. Define x = |fc|(r^ -|-r^)"/^ as a measure of the ratio 
of internal to galactic forces and evaluate it along the unperturbed trajectory. Next, 
integrate the actual orbit to determine whether it remains bound to the center. After 
many samplings of k, one calculates fundix) the fraction of bound particles at a fixed 
value of X- 

Repeating this procedure for different initial positions (5-20rc) and power law shapes 
(n = 1-2) allows comparison of the importance of the various inputs to the calculation. 
Figure |15| displays fundix) for different n. For a fixed ra, the geometric orientation of 
the rocket produces an intrinsic spread in outcomes for 0.2 ^ x ~ 0-8- By comparison. 
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the entire range of n corresponds to relatively small shifts in x- This geometric effect 
is large and irreducible compared to the uncertainties in the initial position, power law 
shape and number of orbital periods. 

To help gauge the importance of internal forces on the binding of orbits, it is useful 
to define the critical ratio of internal to galactic forces by fundiXcrit) = 1/2. While there 
is some variation, Xcru ~ 0.3 is a reasonable estimate. In applying this result to the 
loop dynamics evaluate x = ^r/lV^/'l at apocenter and regard the orbit as unbound if 

X ^ Xcrit- 

These results have been used to speed up the large-scale numerical calculations 
determining the probability of capture and the halo profiles in the sections that follow. 
Essentially, the rocket effect is ignored before capture and the retention and lifetime 
criteria at the epoch of interest are imposed to determine the bound population of living 
loops. Since the retention criterion is generally stricter than the capture criterion little 
error is made and this allows a single simulation to be used for multiple values of Gfi. 




Figure 15: The fraction of radial orbits that remain bound as a function of x foi' power 
law exponents n = 1, 5/4, 3/2, 7/4 and 2 (red, orange, yellow, green and blue, respectively). 
All orbits began at an apocenter of 20 core radii; 100 random orientations for the internal 
acceleration were drawn; for each choice the orbit was integrated many characteristic orbital 
times; the curves summarize the bound fraction. 



6. Network Evolution 

The halo will contain loops born with a variety of positions, velocities, times of formation 
and lengths. This section discusses the birth rate density for loops generated by a scaling 
network. Succeeding sections discuss the average probability of capture for loops created 
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at a given epoch and the net effect when the birth rate density is integrated over all length 
and time scales. 



6.1 Status of Large and Small Loops 



Early simulations of Nambu Goto strings |104| , |105| , ^successfully tackled the 



large scale properties of the network, in particular, the relation of horizon to correlation 
length, characteristic spacing and persistence length. They validated the original idea 



that the network would evolve to reach a scaling solution in simulations ||106|] . The 
attraction of arbitrary initial conditions to a scaling solution insures that the energy 
density in long, horizon-crossing strings never comes to dominate the total energy density 
of the universe. This is a necessary but not sufficient condition to avoid overdosing the 
universe since loops formed during the radiation era must decay by some mechanism to 
avoid a monopole-like problem. For strings that couple only to the gravitational field 
the decay must involve gravitational wave emission. 

Loops are created when the horizon crossing strings are chopped up. It was orig- 
inally thought that loops would form by intercommutations of the long strings at the 
characteristic scale of the horizon (the Kibble one-scale model). Since the early studies 
found clear signs of a gas of small loops and dense kink-filled string segments at the 
simulations' resolution limit it became apparent that some basic understanding had yet 
to be achieved. In a perfect scaling solution all the properties of the network, not just 
those close to that of the horizon, should scale. Characterizing the small scale structure 
and searching for evidence of its scaling has been an ongoing effort. 

The significance of the small scale structure has now become more apparent. It 
turns out that rather than being a detail to be disentangled from resolution issues, the 
small scale structure is inextricably intertwined with the evolution of the network on all 
scales less than the horizon. Small scale structure governs many potentially observable 
features of the network and ultimately has a considerable impact on the expected halo 
clustering properties of loops. This section briefly reviews areas of significance for the 
problem of loop clustering. 

Current numerical simulations of NO strings generate a continuum distribution of 
sub-horizon scale loops including some loops that are much bigger than the resolution 



limit [^, ^ |107|| . Although non-trivial differences exist between the most recent simula- 



tions a point of common agreement is that the abundance of loops near the horizon scale 
is greater than found in the earlier, lower resolution work. Roughly ~ 10% of the length 
of long strings ends up in loops with length 10~^ < l/t < 10^^. For the purposes of this 
paper, these are all "large" loops because they imply Ht » 1 for Gji/c^ « 10^^. A 
second area of mutual agreement is that the population of these large loops is judged to 
be scaling^ and is not a transient artifact. The presence of a population of large loops 



''The extent of scaling and the criteria for scahng differs. See "Note added" in Ref. |^ for a 
comparison. 
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hearkens back to the original expectation that the horizon scale would determine the 
properties of newly formed loops. 

The traditional interpretation of the large loop part of the distribution is that it forms 
via a sequence of multiple encounters in which a horizon-scale loop is cut into smaller and 
smaller progeny. Some of these encounters are with the long, horizon crossing component, 
some with other loops and some are self-intersections in which a loop oscillates so that 
individual parts collide. Ref indicates that self-intersection of larger loops is the 
dominant overall loop production mechanism. An assessment of the clustering of such 
loops is carried out in this paper. 



Both new and old simulations also contain large numbers of small loops. Early 
calculations suggested ||108| , |109|| and recent studies have demonstrated that tiny loops 
are produced in great abundance by the interaction of small scale structure on string 
segments as those segments first form a large scale cusp I^S], The small scale structure 
today owes its existence to non-linear interactions that the string experienced in the past 
at the time it first entered the horizon. Ref. ||3^ quantified the structure in terms of a 
correlator oc {l/tY^ for size scale / at time t and showed that x is completely determined 
by mean network properties (rate of expansion of the universe and rms velocity of the 
strings). The slope of the spatial correlation function of simulations agrees well with the 
theoretically-determined x over the expected range. 



Most loops (measured in terms of number or total length) are created at small 
physical scale with a cutoff set by gravitational damping and theoretically derived to be 
a ~ {Gfi/c^Y'^^^ where x = 0.1 or 0.25 for radiation or matter respectively [^1^, |111| , p9| . 
Numerical simulations do not include gravitational wave damping so it isn't possible to 
verify directly the cutoff prediction. However, the slope of the loop size distribution 
depends upon x, is insensitive above the cutoff to the actual value of the cutoff and 
agrees well with that found in simulation p6|. 



It appears that roughly 90% of the horizon crossing string goes into tiny loops 
||112| , |113|| . There are many factors which will end up influencing the clustering of small 
loops. While there are more small loops than large ones, small loops evaporate more 
quickly than large ones. Loops at the gravitational wave cutoff are ultra- relativistic |^ . 
Time-dilation and cosmic expansion effects (eq. ^.10| ) in the ultra-relativistic limit imply 
Ht ~ {acyGfif^^^^"^ ~ (G^/c2)^^/^^+''\ In the matter era, u = 2/3, x = 0.25 so 
Hr ~ {Gfi/cy-^; in the radiation era, u = 1/2, x = 0.1 and Hr ~ {Gfi/c^f-^^. All 
these results suggest that loops at the cutoff may be unable to cluster but the small 
power of Gfi/c^ means Ht is never very small and the significance of loops larger than 
the cutoff (which live longer) is murky. The multiplicity of factors at play suggests that 
a detailed calculation of clustering should be carried out. 
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6.2 Birth Rate Density 

The long horizon crossing strings are chopped into loops at a rate 

'dp\ 2pc 

oo^loops 

where Poo = /w/ ilst)"^ is the energy density in long strings with typical separation and 
mean square velocity -u^ ||114| . 



) =^(1--(1+C)) (6.1) 

/ OO^IOODS 



As long as a scaling solution is achieved when horizon oc t (true for power law 
expansion in radiation and matter dominated eras but not applicable to the recent A- 
dominated phase) then the average birth rate of loops of physical size / per physical 
volume V at time t 

dN _ fix) 
dldtdV ^ ■ ^ 

where x = l/t = a/{Ht) for some function f{x). 

The loop formation process involves interactions within the network and, simultane- 
ously, stretching of the horizon crossing strings and expansion of the universe. At least a 
few expansion times (~ ^/ H) are needed for intercommutations to transform long string 
segments into sub-horizon loops. Once sub horizon loops are formed the probability for 
loop-loop interactions decreases and the loop achieves a fixed physical size in a few more 
expansion times. The assumption in this paper is that during each infinitesimal time 
interval {t,t + dt) the network produces the loops implied by dN / dldtdV . The loops have 
physical scale that changes only due to gravitational radiation; they suffer no further 
intercommutation. This prescription provides the distribution of initial conditions for 
the clustering calculation. 

In principle, the full joint distribution of dN / dldtdV , loop center of mass momentum 
and rocket direction in the center of mass frame is needed to realize the initial conditions 
for the dynamical motion of a population of loops in the background FRW cosmology. 
In lieu of a detailed description, assume a factorized form for the joint distribution 

dN _ dN dP dP 
dldtdV dvdhr ~ dldtdV dv dClr 

where v is the loop center of mass momentum and fir the direction of the rocket. Here, 
dP/d{..) means the differential probability for the initial condition of (..) with unit nor- 
malization. In practice, the correlations between the magnitude of the momentum v and 
l/t is retained but all correlations between the direction of the momentum, the direction 
of the impulse and l/t are discarded. 

In the homogeneous limit number A^ and length L of loops created in comoving 
volume V = Vq/oq (i.e. the comoving volume implied by a physical volume Vq today; 
ao = a(to)) are 

'n\ Vo fdy' fa'Y f ^ ,( fix') 



L I tlJ y'^ \aoJ J "^"^ [fix')x' ' ^^'^^ 
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where y = t/to, x = l/t and a' = a{t'). By the current time some loops have evaporated 
and all have shorter lengths. The distribution of loops with length is 



di 



where L = Llif , I' , Qf' , t' ; to] is length at to in terms of the initial loop parameters and, 
similarly, 9a = QaIv", I', flr',t'] to] is or 1 depending upon whether the loop has reached 
the end of its life or not. The dependence of L and on dynamical variables may be 
traced to relativistic effects that link FRW and loop center of mass frames. Quantities 
and La are defined by integrals over dN^/dl. 

When the loop's center of mass motion is only mildly relativistic then the loop 
lives until tujf, = t' + V /{GiiVe)- Hence 6^ = 1 for t < tuff, and the loop length 
is L = /' — TEG^{t — t') independent of the loop dynamics. To summarize: for non- 
relativistic kinematics and ignoring the dynamical influence of the rocket the forms for 
L and 9a simplify: L = L[l',t';to] and 9a = Qa[1' ,t';to] giving 



dNA Vo fdy' 



dl tl J y 



(6.6) 



The integration over x' = I'/t' follows directly since occurrences of /' are now easily 
rewritten in terms of /, to and t'. 

6.3 Fragmentation and Cusp-Mediated Loop Formation 

Until this point, the network evolution has been treated in an agnostic fashion with 
respect to the scaling solutions for matter and radiation eras. For the purposes of pre- 
senting numerical results, the focus tightens to matter era network evolution. This choice 
avoids any potential inconsistency with respect to the growing galactic perturbation for 
which matter era dynamics are most appropriate. It also simplifies and streamlines the 
presentation. However, most loops in the halo today were born before equipartition and 
the absolute numbers of such loops depend upon the radiation era expansion dynamics 
0. This paper concentrates on enhancement, i.e. the ratio of clustered to homogeneously 
distributed loops which is expected to be less sensitive to expansion dynamics. 

Large loops are created by hierarchical fragmentation; small loops by cusps interact- 
ing with pre-existing small scale structure. The processes are both active at the same 
time. The description for each mechanism begins with the power law form 

, , I Ax~^ if Xl < X < Xn , , 

fix) = { \ ^ (6.7) 

[ else 

which depends upon A, /3, xl and xu- Individual /'s for individual mechanisms are 
weighted by the fraction 6 of the total energy loss rate by horizon crossing string 
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{dp/ dt) oo-,ioops that ends up in loops created by a given mechanism. Energy balance 
gives 

(6.8) 



7s y^u ^1 

In the non-relativistic matter era, u = 2/3, 7^ = 10^^ and v'^ = 0.35. 

Refs. ^ explored the cusp-mediated mechanism and combined theoretical 

arguments and simulation-derived average quantities to infer in the matter era 

X = 0.25 (6.9) 
f3 = 3-2x (6.10) 

= ^ = 20. (G/x)'+'^ (6.11) 

For the cusp-mediated mechanism / is strongly tilted to small scales; the upper cutoff 
xu has little effect. Ref. |3^ derived the means square velocity distribution for newly 
formed loops 

(/3^) = 1- ... , '1 , fiV (6.12) 



(2x + l)(2x + 2) Vt. 
where B = 0.61. 

In this paper, the theoretically derived / is adopted to describe cusp-mediated loop 
production. The cutoff xu is adjusted freely. Ref. noted that a puzzling discrepancy 
exists between the above expression for and previous, simulation-derived dispersions 
[rT3] with the theoretical result being larger. Since loops with higher velocity must expe- 



rience more damping before they are able to cluster adopting the theoretical expression 
gives the "worst case" scenario for small loop clustering. The implications of a reduced 
will be investigated as well. 

Ideally a similar approach for the hierarchical formation mechanism should be fol- 
lowed. Fits for / for large loops in an expanding cosmology are not generally available 
and a systematic comparison of existing simulation results is lacking. One network sim- 



ulation ||107|| gave a scaling, unimodal distribution at x ~ 0.1 for / but another |^ 
lacked the peak and generated an approximate power law form for / at smaller x. 

A potential practical complication is the extent the cusp-mediated contribution in- 
terferes in fits to / designed to characterize the fragmentation mechanism. Figure 3 in 

showed that cusp-mediated loop production traces x^/(x) from the simulation ||107 



over the range x < 10"^. This subsumes a significant part of the loop range termed 
"large" here. 

In sum, the form for / for the fragmentation mechanism for large loops has not yet 
been sufficiently well-characterized to yield specific values for 6, P, xl and xu-^ The 



*Nor is it clear that a truncated power law form will ultimately be sufficient to trace /. The segment 
of curve 10^^ < x < 10^^ in Figure 3 ([pq]) may not be well fit with a simple power law. 
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values of S and Xu are the most important and the most uncertain input for many 
purposes. 

In this paper S, P, xl and xu are regarded as free parameters for assessing the impact 
on loop clustering. Simulation-based results for (/?^) as a function of l/t ||115|| were fit 



using the same form as eq. |6.12| in a purely empirical manner. In the matter era the 
results are B = 1.46 and x = 0.114. This fit for large loops formed by fragmentation 
will also be employed as an alternative to the theoretically derived (/?^) for small loop 
motions. 

For all mechanisms, the loop momentum distribution (v) is assumed to be thermal 
in the FRW frame, i.e. 

dP ^ dPdP_ 
dv dv dQjj 
dP _ w^e"^^ 

~ f dv'v'^e-^^' ^^'^ ^ 

E{v) = VTTv^ (6.15) 

^ = — (6 16) 

dQi Att 

where k = k,{x) is set by requiring agreement with the fits to the dispersion 

The direction of rocket impulse is assumed to be isotropic in the loop center of mass 
frame. 



7. Loop Clustering 



This section describes the halo profile formed by loops which are born at a single epoch 
and with a fixed l/t. 



7.1 Probability of Capture 

Consider the probability that a single loop formed at time ti in a large but arbitrary 
comoving volume V ends up today bound to the galaxy with physical semi-major axis 
r. Let AA^ be the number of loops formed in infinitesimal time interval t' to t' + dt' for 
t' = U. 

First, write out the formal probability that the loop has not evaporated and is bound 
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with position x, momentum v and length / today 

dxdvdl AN 

1 r , ,„ dPdP dP 



(7.1) 



r dPdP dP 

/ dx'd^fdn,>——-—5'{x-X)6'{v-V)5{l - L)eAB (7.2) 
J dl dv dilfi 

dV ~ Jfdl~ t'F ^'^■^^ 
F = ! f{x)dx (7.4) 



V 

dP f f 



The initial variables are position x' (assumed homogeneous), velocity v', length Here 
X = X{x', v', I', Qr', t'; to) is the formal time-dependent solution for position, likewise for 
V and L. The function 9ab is 1 if the loop is alive (has not yet evaporated) and bound 
to the perturbation and otherwise. The probability that the loop has initial length I' is 
dP/dV and similarly for the other initial variables. Here, F is a normalization constant. 

Evaluate this integral by Monte-Carlo methods for t' = ti\ first, sample Z', x' and v' 
and then use direct numerical integration to evaluate the final positions and momenta. 
Finally, marginalize the multi-dimensional distribution and focus solely on semi-major 
axes r and current loop length /: 

^^^^ ^ dxdv^^Ar - R) (7.5) 



drdl J dxdvdl 

^^^^ ^ dl'^ (T.6) 



dr J drdl 

where R = i?(x, t7, to) is the formal expression for the semi-major axis in terms of the 
current phase space coordinates. A more detailed description is given in Appendix |^. 

The total probability the loop is bound today is proportional to l/V. As a basis for 
comparison, consider the case of a cold dark matter particle: it is judged to be bound 
if it lies within today's comoving turn- around volume Vta- Scale the differential and 
cumulative forms by the same factor: 

dQAB _ V dPAB /„ „N 



QM<r) = l^^^dr i7.i 



'0 

The combination nVtaQAB{< is the expected number of bound objects today within 
distance r for a mean homogeneous density n and turn-around volume Vta- By construc- 
tion, cold dark matter has Qab{< fta) = A^(A = 1) = 97r^/16 = 5.55. The fact that 
Qab{< rta) > 1 shows that the perturbation attracts distant particles so that the total 
within today's comoving turn-around volume is larger than the number in an equivalent 
volume far from the perturbation center. 
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The mean interior density is 

Qab{< r) 

riAB = (7.9) 

also a useful measure of the clustering and highlights the part of the distribution near 
the center. 

7.2 Halo Profiles 

Integrating from t' — ti to the present gives the expected differential number of loops in 
comoving volume V — Vq/qI today 

d^AB _ f dt' fa'V dPAB . . 

dxdvdl ° J {t'y \ao) dxdvdl ^ ■ ^ 

-V,^fJ (^(^) (7.11) 
dNAB p f dt^ ( af_y dQAB 

drdl " 7 (tr [aoj drdl ^ 

The differential probabihties include a factor 1/Ft'. 

Objects must not have fully evaporated (©^ = 1) and must be bound (©^ = 1) 
to the perturbation to contribute to these distributions. The length is integrated as a 
variable so its straightforward to evaluate Qa- In practice, deciding whether a loop is 
bound amounts to checking whether the orbit has experienced multiple passages through 
the perturbation center (capture). If captured, the integration is suspended but it must 
still be determined when the rocket effect detaches it. If the loop is captured and not 
detached at the epoch of interest it is called bound. 

7.3 Truncation by Rocket 

A loop remains bound with approximately fixed semi-major axis until the internal ac- 
celeration exceeds that due to the gravitational potential. The condition for escape is 
X > Xcrit at which point the loop leaves quickly on an orbital timescale. 

When the rocket is ignored, the geodesic trajectory is independent of / and fi i.e. 
just the motion of a test object.^ As the loop shrinks, the rocket acceleration grows 
monotonically oc /x// while the acceleration at apocenter for a captured loop is constant. 
Consequently, eventually x > Xcrit- For loops with non-relativistic velocities 

X - ^ (7.13) 
a 

a. = r,^ = pl-^ (7.14) 

^ ^ ~:?r~ = —GpH{to)rta^ (7.15) 



r 



2 

ap 



®The probability distributions dPAs/drdl and dQAs/drdl depend upon I not just because of / but 
also because of the correlation between the initial velocity and loop length. 
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where r^p is the physical apocenter, the loop lives until tuje = t + l/{GfiTE) and A = 
fap/fta- In practice, t, Tap and / are determined at the instant the loop is judged as 
bound so that the non-relativistic approximation is good. 

The acceleration a oc M.x/^^ is monotonically decreasing with A or Tap- For bound 
loops with given tufe (i-e. time of formation, size and tension) there is a single value for 
Tap today which satisfies x = Xcrit- Denote the solution rap^crit' J' , f^)- 

Assume an ejected loop instantaneously leaves the galaxy. Loops with r^p > Vap^cr 
have been lost; loops with r^p < Tapper are still in the halo. To summarize: 0^ contains 
a factor 9 (r^p < r^p^cr) which accounts for loop loss by the rocket. 



7.4 Results 

The cumulative distribution Qab{< was calculated for specific choices of fi, formation 
time t' and length I' = ac/H' and results are displayed in figure 16-18. Recall that 



Qab{< r) measures the expected number of bound objects within distance r for a mean 
homogeneous density equal to one object per turn-around volume Vta- 

The figure's black line is Q for coUisionless cold dark matter as calculated in the self- 
similar radial infall model. It serves as a standard of comparison for the more complicated 
infall scenario in which loops need to slow down to be captured and eventually are ejected 
by the rocket effect. The colored lines give Q's for loops with different formation epochs. 
All have string tension Gfi = 10~^^ and length I' = a/H' for fixed a = 10"^. Loops are 
a fixed fraction of the horizon size at birth; older loops are smaller. The initial velocities 
were drawn from a fixed thermal distribution with = 0.2 (this is arbitrary and not 
directly tied to any of the string network estimates). 

Early formation gives a profile that closely matches the cold dark matter one at 
small r/rta (leftmost lines). Such loops have had plenty of opportunity to damp by 
cosmic drag so they cluster just like cold dark matter. Note the empty circle at which 
some curves end. The profile is truncated at larger radii because the rocket effect is able 
to unbind orbits at apocenter. To summarize: for fixed a, the oldest loops are smallest, 
closest to the end of their lives, feel the largest rocket effect and may be retained only 
by the centermost parts of the potential. 

For loops that are not as old, outer regions of r/rta are accessible. Note that many 
curves end near r/rta ~ 0.1 without an empty circle. The endpoint is not a consequence 
of physical ejection but of the minimum time needed to bind an infalling loop to the 
perturbation. Unlike a cold dark matter particle which is known a priori to be bound, a 
loop is judged bound only after it has passed back and forth several times through the 
center. 

With a limited number of particles the Monte-Carlo calculation always has an in- 
nermost radii. However, arbitrarily small velocities are present in the initial conditions, 
so arbitrarily small turn-around radii are possible. 
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log 11= -13, log a= -1.00 




Figure 16: The cumulative distribution Q as a function of r/rta- The black line is Q for 
collisionless cold dark matter according to the radial infall model. It includes material infalling 
for the first time as well as regions of multi-streamed flow {r/rta < 0.36). The colored lines 
show Q for loops formed at different epochs which are known to be bound to the perturbation. 
All have G/i = 10^^^ and a = 10^^ and formation lengths = a/H'. The time of formation 
varies t'/^o = 10^^'^ — 10~^; the color key is log^^o^V^o written on the right. Note, for example, 
that the orange line closest to the center is formation at the earliest epoch while the red line 
at larger r/rta is the most recent. These indicate the degree of clustering in loops compared 
to the cold dark matter case. (1) Old loops closely track the cold dark matter. (2) Empty 
circles indicate the truncation of the halo's loop profile because r^p > rapier, i-e. the rocket 
effect strips the loops further out. (3) Q is not plotted beyond r/rta ~ 0.1 - not enough time 
has passed to satisfy the criteria that the loop is bound (that it pass through the perturbation 
center several times). (4) Profiles which terminate at small r/rta do so only because of the 
limited number of particles used in the calculation or the limited extent of the figure; actual 
profiles extend to the center. 



Compared to figure |T^, figure has smaller a = 10^^ while figure has larger 
Gfi = 10~^^. To the extent that the rocket effect is ignored prior to ejection the Q's are 
the same (all were constructed from the same data). The only impact of these changes 
is to shorten the loop lifetime, so removing some of the lines, and to decrease Tapper, 
truncating the profile at a smaller radius. 

Figure |19| presents the mean interior density n for the string loops compared to that 
of the cold dark matter. Evidently, old loops with small /i have sufficient time to cluster 
strongly. 
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8. Current Halo Profile 



The actual halo profile is more complicated than the examples constructed in the previous 
section because it involves loops of many sizes created over a range of epochs. Results 
for two specific formation scenarios are illustrated in this section. 

8.1 Measures of the Loop Distribution 

The cumulative distribution Nab{< r) = J dl drdNAB / dldr and the length-weighted 
cumulative distribution Lab{< = / ^dl drdNAs/ dldr provide summary information 
about the number and total length of all loops bound to the galaxy. These quantities are 
normalized with respect to Na and La, the equivalent quantities expected to be present 
in homogeneous space (eq. [6.51) . A measure of the cumulative number of loops bound 
to the galaxy compared to the total within the turn-around volume is 

If all loops behaved dynamically like cold dark matter particles then one would expect 
Q(< r) to resemble Qab{.< f) (eq. p^.S] ). 

In a similar manner, start with the average number density of alive, bound loops 
within radius r 

_ , . Nab{< r) 

^^^(^) = y(r) ^^^^ 

and the average number density of alive loops in homogeneous space fiA = Na{< r)/V{r). 
A measure of the average overdensity of loops bound to the galaxy is 

m< r) = (8.3) 
riA 

To the extent that loops behave like cold dark matter then one expects A/'(< r) to vary 
like Qab{< r){rta/rf. 

Let X be the average over the length distribution of some function; denote by X^*^ 
the length- weighted moment of the function times T for i = 1, 2,... The first moment 
corresponds to total length or energy. 

The substitutions Nab -^ab ~ ^ab and A^^ N^^^ = La leads to cumulative 
and density with respect to energy instead of numbers. In an obvious fashion, let Q{< 
r) Q(i)(< r) = If^, UABir) ^ n%, ua ^ n«, and N{< r) ^ m^\< r). 

8.2 Point: Large Loops from Fragmentation Model 

A model in which all the long string length goes into sub-horizon-scale loops will be 
discussed first. The model parameters are 5frag = 1, = 10~^, au = 10"^, fi = 10~^^- 
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10 ^ with /(x) <x. X ^ . Let g{x,t) be integrated over the birth rate density: 



dl I dt- 



-a{tyg{x,t) 



oc J dxx'^ J dtt^''"^g{x,t) 



U) 
^.5) 



For g = 1 the integral I[g] is proportional to total number of loops born. For matter 
(radiation) era u = 2/3 (1/2) the result varies like 1/t (1/t^/^) and the number of loops 
is dominated by the earliest epochs. For g = I = xt, I[g] is proportional to total length 
of loops born and the result varies like logt (l/t^^^). For P < 1 large x loops dominate 
as measured by number and by length. 

The cumulative number Q is shown in figure ^for /i = 10~^^ with and without the 
cutoff imposed by the gravitational wave recoil. Qab for cold dark matter (the black 
line) provides a point of reference. All three cumulative distributions are close at small 
radii (r/r^^ < 10~^'^ or approximately < 3 kpc). The rocket is effective at depleting the 
old and hence small loops that would otherwise be present at large radii. Since loop 
numbers are dominated by formation at the earliest epochs, the characteristic signature 
of the recoil effect is the depletion of large numbers of loops at large radii. The inner 
regions are not immune to the depletion but it is less severe. 




Figure 20: 



Clustering for Gfi 

-3 in-l> 



10~^^ by number of loops for / constant (/3 = 0), loop size 
I = a/H, {aL,0!u) = (10^'^, 10^^), matter era expansion dynamics (a oc with z/ = 2/3). The 
black line is log Qab (< ''') for cold dark matter, the blue lines are log Q for loops with (lower) 
and without (upper) the gravitational wave recoil (the rocket effect). 

By contrast Q^^^ weights loops by today's length. This distribution samples a broader 
range of times. The three cumulative distributions (with and without rocket and cold 
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dark matter) now have a very different set of relationships. There are two quahtative 
observations: recoil makes little apparent difference and the amplitude of the cumulative 
loop distribution lies below the cold dark matter one. 



I r 
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Figure 21: Clustering by energy of loops (same parameters as figure ^0|). The black line 
is logQAB{< t) for cold dark matter, the blue lines are log Q*^^) for loops with (lower) and 
without (upper) the gravitational wave recoil (the rocket effect). 

Recoil is not dramatic in Qp''> because only a small contribution to total length is 
made by loops near the end of their lifetimes. This can be understood by revisiting figure 
H which displays the bounds on formation time and string tension. At fixed tension the 
logarithmic interval between the lifetime (red diagonal) line and the current epoch is 
proportional to total loop length created and present in homogeneous space. Only part 
binds to the galaxy, the interval between the lifetime line and the capture time (turquoise 
horizontal) line. The rocket cuts out the space below the acceleration condition (green 
diagonal). For small G/i many decades lie above the acceleration line and below the 
capture time so it is difficult to see the effect of the rocket on the length-weighted loop 
distribution. 

Insensitivity of Q*-^-* is not a prerequisite for clustering but an interesting consequence 
of the weighting and scale factor. One expects the rocket to have a more visible influence 
on Q'^^) in the radiation era when the distribution is oc not logt. 

The second observation, that the amplitude of the cumulative distribution lies below 
the cold dark matter comparison, is related to the necessity of slowing down enough for 
capture. Again, referring to figure H, at fixed G/x loops born between the capture line 
and the current epoch are present in homogeneous space but cannot bind to the galaxy 
today. The effect lowers the amplitude relative to the cold dark matter scenario where 
there is no such constraint. At first glance, the figure would suggest that the amplitude 
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be diminished by ~ 2 (i.e. at Gfi = 10~^^ about half the decades he above the capture 
hne and half below). The figure illustrates Vi = 0.1 whereas the simulation velocities are 
larger (rms dispersions are 0.5 — 0.9). The excluded region in the figure increases and 
accounts for the factor of 3 — 4 diminution in amplitude observed in the simulation. 

As the tension increases the natural expectation is that the curves should fall away 
from the cold dark matter analog since Ht diminishes and less damping will be possible. 
In figure this corresponds to trying to move to the right and the available phase space 



shrinks. The cumulative number Q are shown in figure for range of string tensions 
10~^^-10^^°. The lowermost profiles have tensions Gfi = 10^^^ and 10^^^ and fulfill this 
expectation. However, a striking feature is that all curves with ^ 10^^^ are bunched 
together. They track the cold dark matter profile at small radii and are stripped beyond 
a characteristic radius. Without the rocket effect the same subset of curves closely tracks 
the cold dark matter profile to about r/rta ~ 0.1 (at which point deciding whether or 
not a loop is bound is problematic). 

Since most loops are created at early times, in the absence of dynamics most have 
ages ~ to- The gravitational wave acceleration of such loops is qr ~ Tpc/TEto- Setting 
this equal to Xcrii|V0| gives the characteristic radial scale XcritVc^ Eto / {cT p) ~ 2.3 kpc 
or logr/rta ~ —2.7, approximately what is observed in the simulations and independent 
of Gfi/c^. 

1 r 



logior/^ta 




Figure 22: String tensions G/i = 10"^'' to 10~^^ in powers of ten, bottom to top with all other 
parameters the same as figure |20| ). The black fine is logQAB{< f) for cold dark matter, tfie 
colored fines are log Q (all with rocket effect). For r/rta < 10~^'^ the group with G/.i < 10~^^ 
are numerically indistinguishable given the finite Monte-Carlo sample. 



Profiles with loops weighted by length are shown in ^ Because the loop distribution 
is logarithmic, the measures are not dominated by the loops formed at the earliest times 
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nor is the rocket effect a dominant influence on tlie sliape. As tension decreases, tlie 
cumulative approaches the cold dark matter limit albeit slowly. 




Figure 23: Clustering by energy of loops with string tensions Gfi = 10~ to 10~ in powers 
of ten, bottom to top (colored lines) and all other parameters the same as figure |20| ). The 
black line is log (5a_b(< ''') for cold dark matter, the colored lines are log Q^^-* for loops (all with 
rocket effect). 

The average interior density within a radius is displayed in figures (number) and 
^ (length). They illustrate that the loop clustering follows that of the dark matter. 

Figure ^ shows that varying /3, the slope of the loop distribution function, has little 
effect. The significance of this observation is not that the change in j3 is ignorable but 
that it is subsumed by scaling with respect to the homogeneous results. 

All the calculations in this section will be altered by treating the loop production 
during the radiation era with the correct scale factor a{t) oc t^/^. The weighting of the 
birth rate in the comoving volume will shift from a logarithmic distribution to one that 
oc t"^/^. New results will be calculated in the future but one can infer that this new Q*-^-* 
will have properties intermediate between Q (dominated by early time production) and 
qW (logarithmic) calculated thus far. 
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Figure 24: Mean interior number density of loops log]^o-^(< (relative to unclustered loop 
population, including rocket effect). String tensions Gfi = 10^^'^ to 10"^^ in powers of ten, 
bottom to top (colored lines) and all other parameters the same as figure 20). The black line is 
logQAB{< r){rta/r)^ for cold dark matter. The infall model has turn-around radius rta = 1.1 
Mpc. 




Figure 25: Mean interior energy density of loops log;^Q AA*^"^) (< r) (relative to unclustered 
loop population, including rocket effect). String tensions Gfi = 10^^^ to 10^^^ in powers of 
ten, bottom to top (colored lines) and all other parameters the same as figure [20|) . The black 
line is logQABi< f){fta/fY cold dark matter. The infall model has turn-around radius 
rta = 1-1 Mpc. 
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Figure 26: Cumulative distributions for /3=0 (blue), 1 (green) and 1.63 (red) for for Gfi = 
10"^^ for number (left) and energy (right) of loops. The remaining parameters are the same as 
figure pO|. 
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8.3 Counterpoint: Small Loops by Cusp Formation 

According to current estimates 80-90% of tlie liorizon-crossing string ends up forming 
small loops. They are unlikely to cluster well because they are born moving fast and 
they have intrinsically small lengths. However, they have a broad range in x = l/t and 
the relativistic kinematic effects are non-trivial (time dilation, energy shifts). So it is 
of interest to investigate the degree of clustering within the galaxy. For a simplified, 
canonical treatment assume 100% efficiency (5 = 1) for chopping long strings into loops 
characterized by f{x) oc a;~^ over the range x = {xqw, !)• Here, the power law index 
j3 = 3 — 2x and x = 0.25 in the matter era. The loops are weighted in number and mass 
to the low cutoff at the gravitational wave damping scale xgw = 20. (G/i)^^^^. The 
velocity dispersion of newly formed loops is given by the ultra-relativistic theoretical 
result eq. |6.12| . 

Results for G/i = 10~^^ are presented in figure for the cumulative number dis- 
tribution Q and the cold dark matter behavior Qab- It is apparent that only a very 
small fraction of the loops generated by cusp formation are able to live long enough 
to slow down and bind to the galaxy. Recall that homogeneous space has Q = 1; the 
majority of small living loops in proximity with the galaxy are not bound to it. The 
figure also illustrates the effect of turning off the rocket effect and starting from a less- 
relativistic velocity distribution function. These change Q in the manner expected but 
still leaves the relative numbers small. A similar picture emerges from the cumulative 
length distribution Q^^^ in figure 



The average interior number and energy densities are shown in figures |2^ and |30. 
The canonical model reaches the background number (energy) density at r/rta = 10"^'^ 
(10~^). The smallest r/rta at which the infall model will be a reasonably accurate physical 
description of the galaxy is ~ 10~^. 

A comparison of Q^^^ for tensions G/i = 10~^^-10~^^ shows that increasing /i increases 
the Q^^\ Essentially, the gravitational wave cutoff xqw yU^'^ increases and each horizon 
crossing string makes more large loops which live longer and move more slowly. The 
range of uncertainty is shown by comparing the results for velocities drawn from the less 



relativistic Ref. 115 
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Figure 27: The yellow unadorned line is the cumulative number distribution of loops (log Q) 
for G/i = 10~^^ in cusp-generated loop formation during the matter era {u = 2/3, x = 0.25, 
v'^ = 0.35, 1/7^ = 10, {xl,xh) = (xcwA) with xcw = 20./x^+^^ and with theoretically 
determined ultra-relativistic velocity distribution eq. 6.12| ). The black line is logQA_B(< r) for 
cold dark matter. The lines with dots represent altering the standard model assumptions: sup- 
press rocket effect (blue) , adopt less relativistic initial velocity distribution (red) , and combined 
(green) . 
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Figure 28: Cumulative energy distribution of loops (log Q^^^); identical parameters as figure 
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Figure 30: Average interior energy density of loops (logM^^^); identical parameters as figure 
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Figure 31: The cumulative energy distribution of loops (logQ^-^-*) for Gfi = 10^^^ — IQ-^^ 
(yellow to blue) in cusp-generated loop formation during the matter era (z/ = 2/3, x = 0.25, 
■^oo = 0.35, 1/7^ = 10, {xl,xh) = (xgwA) with xqw = 20. fi^^'^^ and with theoretically 
determined ultra-relativistic velocity distribution eq. 3.12). The black line is \ogQyiBi< f) 
for cold dark matter. The lines with red dots show the effect of altering the standard model 
assumptions to adopt less relativistic initial velocity distribution. All lines include rocket effect. 
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9. Conclusions 



Cosmic superstrings are produced towards the end of brane inflation, a scenario in mod- 
ern superstring theory. As string tension /x decreases and/or formation size a increases 
(where a is roughly the fraction of the horizon), then the timescale for the loop to evap- 
orate by gravitational emission r oc increases. This has important consequences for 
the abundance of strings as potentially observable astrophysical objects[^, ||. Although 
born moving at close to the speed of light, long-lived loops damp by the expansion of 
the universe and fall into potential wells created by cold dark matter and baryons after 
equipartitionH, ^. 

This paper explores the dynamics of loops in the vicinity of the galaxy. The model 
used to represent the growing galaxy is a self-similar, radial infall model for the cold 
dark matter component. The form for the perturbed FRW metric then provides the 
playground for the fully relativistic treatment of string loop dynamics. 

A succession of increasingly complex dynamical investigations addresses the basic 
question "how did the halo get its loops?" These include: capture of fast moving objects 
by the growing perturbation, recoil by radiated gravitational waves, protection of tightly 
bound orbits by adiabatic invariance and a critical transition from confined to free, 
liberated motion. 

The basic picture that emerges is straightforward. Loops damp by cosmological 
expansion coming nearly to rest and some find themselves in the vicinity of a growing 
matter perturbation. These fall into the perturbation and, like cold dark matter, acquire 
radial orbits with scale comparable to the turn-around radius. They linger at roughly 
fixed physical radius while both the galaxy and the turn-around radius grow with age. 
A snapshot of the galaxy would reveal the oldest loops near the center and the newest 
ones on the periphery. The oldest loops are also the smallest loops because the horizon 
size oc t. 

All loops shrink by emission of gravitational wave energy and are subject to sig- 
nificant recoil because the antenna pattern is highly anisotropic. The length decreases 
and the non-gravitational acceleration increases in tandem. The specific assumptions 
about recoil made in this paper are conservative in the sense that they maximize the 
importance of the rocket effect and minimize the possibility of loop clustering. 

For loops of a given size and age the galaxy is stripped from the outside in. Con- 
versely, these loops are retained latest in the central parts of the potential well where 
the binding is greatest. Eventually as the length vanishes and the acceleration diverges 
all will be removed. 

The halo is grown from loops of different sizes accreted over a range of times subject 
to the dynamical capture process. Since small loops move faster than large ones at birth 
two different sized loops born at the same time will generally not accrete at the same 
time. The halo is a mix of loops with a range of sizes and ages at a given apocenter. 



- 55 - 



Two broad classes for loop formation were characterized. Loops directly generated 
by fragmentation ("large" loops) and those mediated by cusp formation ("small" loops). 
It was anticipated and verified that only the large loops cluster about the galaxy to a 
significant extent. 

The bottom line results for large loops are presented in figures (density of the 
number of loops within the galaxy relative to the unclustered, homogeneous value) and 
p5| (density of length or energy of loops). There is a substantial degree of enhancement 
in both number and energy density at a broad range of galactocentric radii that depends 
upon string tension. 

The bottom line results for small loops are presented in figures ^ (number density) 
and |30| (energy density). Little enhancement is observed and that only at radii <1 kpc 
where the radial infall model is not applicable. 

Some essential input physics needs improvement: 

• The loop production function needs to be characterized accurately, especially the 
fraction of horizon-crossing strings that form large loops {5 frag) and the character- 
istic scale of the large loops {xu)- 

• The secular changes to the loop parameters as the loop shrinks need to be calcu- 
lated. This determines the timescale for intrinsic variation of the recoil direction, 
the possibility of loop precession (angular momentum changes not along the an- 
gular momentum direction), and the propensity to convert non-self- intersecting to 
self-intersecting loops. 

These have direct implications for the number and size of the loops that are formed, the 
efficacy of the rocket effect and the loop lifetime. 

From the viewpoint of clustering dynamics these areas need attention: 

• The string network and background cosmology need to be calculated in the con- 
text of ACDM cosmology. All calculations here use the Einstein-de Sitter non- 
relativistic matter model. 

• The galaxy formation model should be made more realistic and consistent with 
structure formation in ACDM cosmology. 

• Dissipative effects may be important in two contexts: dynamical friction can slow 
down loops before clustering begins and can remove energy from loops orbiting 
within the galaxy. 

• Loop-loop interactions may develop at the galactic center as nearly radial orbits 
converge. If such interactions do occur then the small loops formed by intercom- 
mutation may acquire relativistic velocities and be ejected. 
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The first two items should generally improve the treatment of the loop distribution and 
bring it up to date vis-a-vis modern cosmology. The third and fourth items are unique 
to the loop dynamics. They may increase and decrease, respectively, the theoretically 
calculated halo densities. A paper on extending the current results to loops born in the 
radiation era is being prepared. 

Prom the perspective of detectable astrophysical signatures a clustered loop halo is a 
natural source of signals for the following experiments: direct detection of gravitational 
wave emission, pulsar timing variation, and microlensing. Experiments sensitive to the 
local population, especially microlensing of Galactic stars, will see the most significant 
impact of clustering. 
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A. Equations of Motion in Inhomogeneous FRW 

The 4-vclocity and direction of the 4-impulse are parameterized in terms of v, v, n, n 
according to 



V^^{ V^.^l (A.1) 



ATM ^ ( ±^^, ) (A.2) 



= Vl + 2$ (A.3) 

av, - ^1 + 2^ (A.4) 

where c = 1 and hat quantities are unit vectors. The sign of is the sign oi n -i). The 
equations of motion are 



dv (n ■ -0) ar-na^ {v-\/ip) \^T^\^v^ va 

dt + v"^ aa^a^ a a\ 



(A.5) 



dt + aalvT+l^ 

vTTv2((w-v^) V* - -0,?) 

dn {n-'Vip) {n ■ v) nv (n ■ v) a^a^v {n-v)'^nv'^a {n ■ v)'^ nv'^^ /A7^ 

dt ~ aa^av-VlT^ ^ v/TTl^ a{l + v^) al{l + v^) ^'^ 

dn" arO£^v{{h ■ v) n'^ — v^) {h ■ v) v'^{{n ■ v) — v^)d / \ q\ 

'dt ~ nVT+l? ^ a(l + t;2) ^ ^ ' 

{n ■ v) n^a% - (n- V<l>) a^v^ -\- {h ■ v) a^^^i - {n ■ v) Qi|'i/',i^ 

(ji ■ c) i'-{{n ■ r) n' - r')<i> 

dx* a^vv'^ 
dt aa^Vl + f ^ 



+ 



(A.9) 



da,,. rE;ar.^a,/, 

- ^ ' ^ (A.IO) 



dt rpVIT 



Terms have been organized into 3- vector dot products {{d ■ e) = J^j^^ 3 ci*e*), the exphcit 
index i runs 1 — 3 and time derivatives are indicated e = ^. 

at 

B. Relating Center of Mass and FRW Accelerations 

Let be the initial 4- velocity of the string center of mass in the FRW frame. For clarity, 
let V^-^ be the components in the FRW frame ( "x" ) and first define a transformation to 
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an orthonormal frame ("y") by 

dy^ = ^/l^^dx^ (B.l) 

dy' = a^/lT2^dx' (B.2) 

Second, introduce a Lorentz boost to the rest frame ( "z" ; designated A^) i.e. hy{y) — > 
= (1,0,0,0). Explicitly 

dy^ 



The impulse in the string center of mass frame is a"^^ = (0, a^^-j) = a^(0, n^^-j) where n^^^ 
is a unit-vector. The 4-impulse in the FRW frame is 



'1+2$T/ 



(B.7) 



arN"". (B.8) 



where V^^) ■ n^,) = Efc=i,3 and V^^) ■ V^^) = Efc=i,3 This specifies the 



initial A^". 

It remains to determine the scalar Ur- In the center of mass of the loop the rates of 
energy and momentum loss are 

^ = (G/.) (B.9) 

ar^TpiGij)-^ (B.IO) 

t{z) 

bmce, dz^ = ^ the length may be straightforwardly expressed in terms of the FRW 



time 



f dx^ 



"init (a;) 

Loops at highly relativistic speeds in the FRW frame have extended lifetimes because 
their center of mass clocks advance more slowly. It is convenient to express the evolution 
of ttr directly in the FRW frame using the global time coordinate 

d f l\ Vi 



dx^ \ar) VpV^ 



(B.12) 



After a^A"" is initially set the entire calculation can be carried out in the FRW frame. 
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C. Variants of Figure 9 



Figure ^ shows the bounds on the formation time and string tension for a loop with 
initial velocity Vi = 0.1 captured and retained at physical radius 30 kpc. It is essen- 
tially the same as |^ but the black lines are based on v found by numerical calculations 
using eqs. ( p.l| )-( pl4D for aligned and anti-aligned rockets rather than by the analytic 
expressions in § |5.2| . It generally validates the analytic approximations, however, it shows 
the existence of some additional phase space for capture when the relative directions of 
impulse and initial velocity vary. The green line is the retention criterion based on the 
critical acceleration at the current epoch. 

Figure ^ is similar to the one above but applies to capture at a smaller physical 
radius 10 kpc. The formation time must be earlier and the acceleration limit approaches 
the loop lifetime limit. 
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Figure 32: Same as ^ except that v is calculated with the full relativistic equations of motion 
for aligned and anti-aligned rockets to determine capture and the numerical experiments for 
randomly oriented rockets are omitted. 
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Figure 33: Same as ^ except for physical radius 10 kpc. The red lines (damping to the 
characteristic rotation velocity and lifetime) are unchanged. All the other lines move. Capture 
requires earlier formation and retention allows larger G//. 
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D. Monte- Carlo Dynamics 



The probability that the particle is bound with position x, momentum v and length / 
today 

dPbnd _ 1 f .-tt.-^.pdPdP ^ 



dxdvdl V 



f dx'difdl'^^5%x - X)S^{v - V)S{1 - L)9ab {x,v) (D.l) 

/ (aiU tilt 



where X, V and L are the formal time-dependent solutions for initial position x', mo- 
mentum if and length /' at time t'. The integration is over the distribution of position 
and momentum while /' = /j and t' = ti are fixed. 

There are some significant computational simplifications. First, only the properties 
of the bound particles with / > arc of interest. To avoid sampling initial conditions 
which imply 6*^ = and/or loops that have evaporated use the following approximations. 
In flat FRW space, directly integrate the motion and length of a particle with given initial 
momentum from to Iq. The particle must live to the current epoch and be able to 
enter the spherical volume defined by today's turn-around radius or there is no chance 
of it being accreted. This limits the minimum length as well as the maximum distance 
and the range of angles between x and v that need to be sampled. For acceptable I' 
the Monte- Carlo points are weighted by the ratio of the volume actually sampled to the 
total volume. 

For efficiency order the kinematic integrals (1) magnitude of v' , (2) magnitude of x', 
(3) direction x' and (4) direction p' . First, sample v' according to dP/dv' (weight is 1). 
Second, calculate Xmaxi the maximum initial displacement from the perturbation center 
that a particle can have and still reach the turn-around volume by today, and sample 
x' from the volume associated with Xmax (weight is A.t:xI^^^/2)V). Third, sample x' in 
the full sphere (weight is 1). Fourth, sample v' from the maximum angular range that 
allows motion from x' to reach the turn-around volume by today (weight is ratio of the 
angular extent sampled to the full 47r extent). 

To generate the final distribution, weight each bound particle according to the sam- 
pling above. In the end, ignore the orbital phase of the particle, i.e. marginalize the 
two-dimensional distribution to give the distribution of semi-major axes 

^ diiv^SS^Sir - R) (D.2) 



drdl J dxdvdl 

where R = R{x,v,tQ) is the formal expression for the semi- major axis in terms of the 
current phase space coordinates. 

The second simplification is to recognize that once a particle is bound to the per- 
turbation the physical dimensions of its orbit are fixed if the rocket effect is ignored. 
Therefore, its unnecessary to integrate all particles from ti to to. Stop once a particle 
has passed back and forth through the perturbation center several times and been cap- 
tured. The error (i.e. the characteristic size of the change in r in the future due to the 
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growth of the perturbation) can be inferred by the variation in the relative heights of the 
peaks in figure ^. To the extent that the orbit is exactly radial, one could extrapolate 
from the last calculated r to the asymptotic one. In this work, the integration after A''^ 
bounces {Nc = 3-8) is halted and the value of r at that time is adopted. A check that 
the final results are insensitive to Nc is made. 

The third simplification (used in §0and^ is to ignore the rocket effect before capture 
and apply the retention and lifetime criteria at the current epoch. Retention is more 
stringent than capture so little error is made. This allows the same simulation to be 
used for different choices of /i. 

A fourth simplification is to work directly with the cumulative distribution rather 
than the differential one. It is ultimately necessary to record only the weight and the 
semi-major axis of the bound particles generated by Monte-Carlo sampling. A sort of 
the semi-major axis of the final sample allows simple construction of the cumulative 
probability distribution from the weights. 
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